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Both  a  coherent  and  incoherent  target  filter  are 
compared  on  an  output  signal-to-noise  power  ratio  basis.  Three 
spread  spectrum  waveforms  are  tested:  intrapulse  biphase 
coded,  interpulse  frequency  hopped,  and  interpulse  frequency 
hopped  with  subpulse  chirp.  In  all  three  cases,  the  incoherent 
target  filter  provided  superior  output  power.  Typical 
incoherent/coherent  signal  power  ratios  were  1-3  dB,  0.75  dB, 
and  0.18  dB  for  the  phase  coded,  FM,  and  FM  with  subpulse 


chirp  respectively. 
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I .  Introduction 
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Problem  Statment 

The  application  of  spread  spectrum  technology  in 
military  radar  systems  may  provide  several  desirable 
characteristics  normally  associated  with  spread  spectrum 
communication  system  techniques.  Spread  spectrum  communication 
systems  are  well  known  and  provide  the  following  benefits  to 
the  communicator:  1)  high  resistance  to  jammers,  2)  low 
probability  of  signal  interception,  3)  simultaneous  use  of 
signal  spectrum  by  multiple  users,  4)  multipath  suppresion, 
and  5)  high  resolution  ranging.  The  feasibility  of  applying 
spread  spectrum  to  radar  has  been  studied  by  Salzman  (Ref'll) 
with  the  emphasis  on  exploiting  the  low  probability  of 
intercept  and  jammer  resistance  characteristics  of  such 
wideband  signal  formats.  The  general  conclusion  of  the  study 
asserts  that  spread  sprectrum  is  applicable  to  airborne  radar 
(utilizing  coherent  pulse  trains)  within  a  decade  or  two.  The 
major  obsticle  is  the  high  sampling  rates  required  to  process 
the  signal.  The  sampling  rate  must, in  general,  be  on  the  order 
of  the  signal  bandwith:  approximately  100-500MHz  for  a  typical 
spread  spectrum  scenario.  The  utilization  of  such  bandwith 
extent  yields  high  range  resolution  properties.  It  is  possible 
that  a  30  meter  target  could  be  resolved  to  1  meter,  resulting 
in  the  appearance  of  30  distinct  returns.  This  target  spatial 
break-up  or  super-resolved  return  can  appear  in  several 
different  range  bins.  Futher,  each  individual  return  may  be 


too  weak  for  detection.  Therefore,  simply  integrating  the 
adjacent  range  bins  (non-coherent  filtering)  would  lead  to 
incomplete  range  collapse  losses.  In  an  effort  to  put  the 
distributed  target  back  together,  a  target  filter  has  been 
suggested  by  Salzman  (Ref'll).  This  is  a  narrowband  ,  matched- 
filter  designed  to  collapse  the  doppler-f iltered  return  signal 
into  a  centriodal  peak  at  the  proper  time-delay  value.  The 
objective  of  this  study  is  to  determine  if  the  target  filter 
is  feasible.  More  specifically,  is  a  narrowband  filter  more 
effective  in  target  range  collapse  than  a  non-coherent  filter. 
The  criterion  for  comparision  is  the  resultant  signal -to-noise 
ratio. 

Scope 

The  study  will  cover  radars  using  several  different 
waveforms  with  extents  on  the  order  of  100  MHz.  Various 
modulations  will  be  computationally  simulated  and  include:  bi¬ 
phase  coding,  pseudo-random  intrapulse  and  interpulse 
frequency  hopping  with  subpulse  chirp  applied  to  the  frequency 
hopping  signals.  Pseudo-random  codes  will  be  generated  to 
control  the  phase  and  the  frequency  variations  and  thereby 
spread  the  spectrum.  Other  waveforms  are  not  considered  for 
reasons  developed  in  chapter  II.  Return  signal  formats  will  be 
constrained  to  be  a  superposition  of  time-  delayed  versions  of 
the  transmitted  signal.  This  is  applicable  to  the  noise-type 
waveforms  simulated  under  the  assumption  of  negligible  doppler 
shift  between  scatterers  (Ref4:52). 


Problem  Solving  Approach 

In  order  to  determine  the  feasibility  of  the  narrowband 
target  filter,  a  computer  simulation  will  be  implemented. 
Signal  parameters  are  specified  by  the  user  and  the  program 
generates  a  pseudo-random  code  sequence  of  length  63  which 
spreads  the  signal  (via  frequency  or  phase  modulation  as 
previously  stated).  Next,  the  user  specifies  one  of  three 
types  of  calculations:  1 )  produce  ambiguity  diagram  plot  and 
signal  spectrum  plot,  2)  produce  received  signal  spectrum  plot 
and  matched-filter  response  plot,  or  3)  produce  received 
signal  spectrum  and  target  filter  response  plots.  The  received 
signal  is  generated  assuming  Swerling  I  or  II  type  target 
cross-section  fluctuations ( RCS ) .  Also,  the  target  filter 
spectral  characteristics  are  assumed  to  be  the  product  of  the 
signal  spectrum  and  the  range  extended  scattering  function. 
According  to  Salzman(Ref 1 1 ) ,  the  range  extended  scattering 
function  has  a  Gaussian  frequency  spectrum.  Various  computer 
runs  will  be  made,  testing  the  waveforms  developed  in  Chapter 
II  and  the  coherent  target  filter  response  compared  to  the 
incoherent  target  filter  response.  The  results  will  be 
complied  and  tabulated  for  each  waveform  and  type  of  RCS 
fluctuation.  The  basic  scenario  consists  of  ten  scatterer 
target  occupying  ten  range  bins. 

Presentation 

This  study  begins  with  an  extensive  review  of  the  study 
done  by  Salzman(Ref 1 1 ) ,  with  emphasis  on  waveform  selection. 
The  various  type  of  modulation  techniques  are  presented  along 


with  fundimental  definitions  and  elaboration  on  the  problem  of 
target  super-resolution.  The  next  chapter,  chapter  III, 
develops  the  necessary  theoretical  tools  to  understand  and 
implement  the  computer  simulation.  In  this  section  are: 
complex  signal  notation,  matched  filter  theroy,  optimum 
receiver  in  white  Gaussian  noise,  received  signal 
characteristics,  development  of  the  radar  ambiguity  function, 
and  a  sampling  theory  review.  Chapter  IV  describes  the 
rudiments  of  the  simulation  software  and  in  this  aspect, 
certain  signal  parameters  such  as  pulse  width,  number  of 
subpulses,  and  time-bandwidth  product  are  presented.  Chapter  V 
presents  the  results  of  the  simulation  with  tables  and  graphs. 
Finally,  ChapterVI  states  the  conclusions  and  recommendations. 


II  Background 


The  successful  operation  of 

an  airborne 

radar 

in 

any 

hostile  electronic 

warfare  scenario 

will  depend 

on: 

1  ) 

the 

radars  ability  to 

operate  without 

being  detected; 

2) 

the 

capability  to  resist  electronic  counter  counter  measures (ECM) 
(Reflljll).  ECM  consists  of  any  electronic  technique  which 
disrupts  the  operation  of  the  radar  system.  Examples  include 
noise  or  spot  jammers,  repeater  transponders,  chaff,  and 
decoys.  The  ability  to  resist  detection  is  generally  termed 
low  probability  of  intercept { LPI ) ,  while  ECM  resistance  is 
termed  electronic  counter-counter  measures ( ECCM ) .  Salzman 
(Ref11:34)  has  concluded  that  coherent  airborne  radar  systems 
can  optimize  both  LPI  and  ECCM  capabilities  via  spread 
spectrum  techniques.  The  report  also  indicates  that  the  most 
significant  advantages  of  a  spread  spectrum  radar  system  can 
best  be  utilized  by  the  fighter/attack  mission(Ref 1 1 : 34 ) . 

The  term  spread  spectrum  is  frequently  used  to  describe 
systems  or  signals  which  possess  one  of  the  following:  1 ) 
signal  modulation  schemes  with  an  RF  bandwidth  considerably 
wider  than  required  for  normal  operations  or  2)  waveforms  with 
time -bandwidth  products  on  the  order  of  100  to  1000.  Spread 
spectrum  communication  most  often  means  the  transmit  bandwidth 
is  several  orders  of  magnitude  greater  than  the  data 
bandwidth.  Spread  spectrum  radar  implies  that  the  transmit 
bandwidth  is  several  orders  of  magnitude  greater  than  the 
reciprocal  of  the  desired  range  resolution.  Regardless  of 


whether  the  spectral  extent  is  much  greater  than  required  for 
either  data  transmission  or  range  resolution,  the  critical 
feature  of  any  spread  spectrum  system  is  the  broad  transmitted 
bandwidth.  These  bandwidths  may  range  from  tens  to  hundreds  of 
mega-hertz.  It  is  useful  to  note  that  the  second  definition 
above  is  also  referred  to  as  pulse-compression.  Historically, 
pulse-compression  was  a  result  of  energy  conservation  efforts 
by  radar  engineers.  For  clarity,  spread  spectrum  will  mean 
that  the  transmit  bandwidth  is  much  larger  than  the  data 
bandwidth(ie  definition  1  above ) (Ref 1 1 : 9-1 0 ) .  These  large 
bandwidths  will  dictate  special  system  hardware  requirements 
and  prudent  radar  waveform  selection. 

The  design  of  a  spread  spectrum  radar  system  begins 
with  the  selection  of  a  signal  waveform  and  it  is  this 
selection  which  the  system  design  must  accomodate.  The 
paramount  criterion  for  selection  of  a  waveform  will  be  the 
maximum  LPI/ECCM  benefit  given  acceptable  ease  of  system 
implementation,  complexity,  and  cost.  The  spread  spectrum 
waveforms  discussed  are  categorized  according  to  their 
respective  modulation  techniques. 

A  RF  carrier  may  be  modulated  via  amplitude  or  angle 
(  -  or  frequency)  and  is  represented  in  comlex  form  as 

s(t)=u(t)expj[2n<f#  ♦f,(t)  )♦  $<t)]  (1) 

where  u(t)  is  the  signal  envelope,  f # i s  fhe  carrier  frequency, 


<*><  t )  the  phase  modulation,  and  yt)  is  the  frequency 

modulation.  The  real  part  of  s(t)  constitutes  the  physical 
signal.  Typical  of  coherent  airborne  radar  are  pulse-doppler 
systems  which  can  be  grouped  as  either  intrapulse ( internal  to 
the  pulse)  or  interpulse( pulse-to-pulse )  modulation.  Further 
modulation  classification  would  be  discrete  verses  continous 
changes  in  phase,  frequency,  or  amplitude.  For  a  series  of 
pulses  the  complex  signal  is 

s<fc>a  H  un(t„ -T)expj[2n(f0t+fri  <t„  -T))+<t»n]  (2) 

where  the  summation  is  not  a  mathematical  one,  but  rather  a 
description  of  a  series  of  pulses.  The  subscript  n  above 
indicates  discrete  coding  by  either  f requency ( fn )  or  phase 
( <t>„  ).  Imbeded  in  un(t)  is  any  continous  modulation  internal 
to  a  given  pulse  and  can  further  be  expressed  as 

un(t)=  An  ( t)exp[  je(  t )  ]  pj[t)  (3) 

where  0(t)  is  continous  angle  modulation,  An  (t)  amplitude 
modulation,  and  p„(t)  is  the  pulse  function  defined  by 

p„(t)=  rect(!l2L(i  )  =  u(t)-u(t-  );interpulse 

rect(ll Zli  ) =  u(t)-u(t-T)  jintrapulse 
T 

where  u(*)  is  the  unit  step  function.  The  difference  between 
interpulse  and  intrapulse  schemes  is  illustrated  in  figure  1. 


Salzman(Ref 1 1 :80 )  has  designated  the  wave forms (modulation 
techniques)  shown  in  table  I  as  candidate  for  the  coherent 
airborne  mission.  In  order  to  select  one  of  the  signal  formats 
of  table  I,  system  implementations  must  be  reviewed  to  insure 
the  feasibilty  of  any  given  selection.  This  has  been  done  by 
Salzman(Ref 1 1 :82-1 09)  with  the  following  modulations  selected 
for  further  consideration:  random  FH,  bi-phase  intrapulse 
coding  with  dwell -to -dwell  PRF  stagger.  The  rudiments  of  the 
above  implementations  will  be  discussed  next.  The  system 
configurations  are  offered  for  completeness  and  should  not  be 
viewed  as  the  only  realization. 

Varying  the  frequency  from  pulse  to  pulse  ,  as  shown  in 
figure  2,  constitutes  interpulse,  random  frequency  hopping. 
The  transmitted  signal  will  be  wideband  but  the  processing  can 
be  either  wideband  or  narrowband.  Wideband  processing  will 
yield  high  range  resolution  and  at  the  same  time  require  the 
frequency  hop  to  be  commensurate  with  the  pulse  bandwidth. 
This  will  help  to  reduce  clutter  and  range  ambiguous  targets. 
A  particular  system  is  shown  in  figure  3.  This  implementation 
has  as  its  main  advantage  an  A/D  converter  which  requires 
sampling  rates  that  are  currently  realizable.  In  other  words, 
the  bandwidth  at  the  output  of  the  phase  detectors  is  on  the 
order  of  the  pulse  envelope  bandwidth;  thus  requiring 
conventional  sampling  rates  of  the  A/D  circuitry.  Phase 
shifting  and  combining  samples  of  the  pulse  train  will  restore 
the  spread  spectrum  resolution.  Physically,  this  is 
accomplished  in  the  range  bins  which  correspond  to  a  given 
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doppler  shift.  The  process  is  repeated  for  each  sample  and 
once  the  range  bins  have  been  processed  for  all  samples,  the 
entire  range  process  is  repeated  for  each  doppler  bin.  Doppler 
bins  are  seperated  in  frequency  by  the  reciprocal  of  the  code 
length  produced  by  the  code  generator. 

Intrapulse  modulation  consists  of  a  pulse  train  where 
each  individual  pulse  is  modulated  discretely.  This  modulation 
can  be  viewed  as  a  series  of  subpulses  whose  phase  or 
frequency  content  is  dictated  by  a  pseudo-random  code.  The 
utilization  of  bi-phase  coding  or  frequency  hopping  for  the 
subpulse  modulation  yields  a  maximum  ECCM  capability  at  a 
minimum  in  target  degradation  and  clutter  rejection. 

Biphase  coding  consists  of  a  transmitted  pulse  with 
subpulse  shifted  fl  radians  in  phase  from  the  center 
frequency.  The  phases  of  each  subpulse  are  determined  by  a 
known  code.  The  frequency  and  temporal  characteristics  are 
illustrated  in  figure  4.  The  generating  code  is  N  bits  long 
and  is  determined  based  on  its  affect  on  the  waveform 
ambiguity  function.  The  ambiguity  function  is  reviewed  in 
Chapter  III.  Each  subpulse  represents  one  code  bit  and 
requires  B=1/T  Hz.  The  critical  feature  is  that  the  bandwidth 
and  hence  the  range  resolution  is  determined  by  the  fine 
structure  of  the  temporal  signal.  Comparatively,  an  uncoded 
pulse  would  require  B/N  Hz,  thus  biphase  coding  results  in 
spreading  the  spectrum  by  a  factor  of  N.  The  corresponding 
range  resolution  is  r=  1/2  cT,  where  c  is  the  speed  of  light, 
r  the  range,  and  T  the  time  extent  of  each  subpulse.  An 


example  of  one  particular  processing  scheme  is  shown  in  figure 
5.  Since  each  transmitted  pulse  contains  the  full  frequency 
spread,  the  entire  bandwidth  must  be  procesed  upon  return  of 
the  signal  in  order  to  insure  complete  correlation.  This  in 
turn  requires  the  A/D  circuitry  to  sample  at  a  rate  as  high  as 
the  signal  bandwidth.  Therefore,  biphase-coded  spread  spectrum 
is  currently  limited  to  approximately  100  MHz  or  less  by  the 
A/D.  The  most  effective  solution,  in  general,  would  be  to 
relieve  the  programable  signal  processor  (PSP)  from  processing 
the  entire  bandwidth.  One  such  method  might  be  a  "target 
filter"  placed  before  the  PSP,  which  would  coherently 
narrowband  the  return  (pre-detection  range  collapse).  The 
target  filter  will  be  the  focus  of  this  paper  and  explained 
more  fully  in  later  chapters. 

Pseudo-random  intrapulse  frequency  hopping  is  similiar 
to  intrapulse  phase  coding  as  previously  discussed.  The  only 
difference  is  that  now  each  subpulse  is  frequency  shifted 
according  to  the  pseudo-random  code.  Figure  6  reveals  the 
waveform  modulation  characteristics.  Figures  7  and  8  depict 
two  processing  techniques,  the  first  is  essentially  a  matched 
filter  receiver  and  the  second  is  termed  "stretch"  processing. 
The  matched  filter  receiver  uses  tapped  delay  lines,  filters, 
phase  shifters,  attenuators,  and  a  coherent  summer  to 
correlate  the  return.  The  stretch  receiver  mixes  the  return 
with  a  local  oscillator  that  is  stepped  in  frequency  according 
to  the  transmitted  code.  The  result  is  that  the  processed 
signal  will  have  a  frequency  shift  that  corresponds  to  the 


target  range.  The  siganl  is  sampled  by  the  A/D,  Fast  Fourier 
transformed  to  form  range  filters,  and  finally  doppler 
processed  and  range  collapsed  prior  to  detection.  The  stretch 
techniques  primary  is  that  it  removes  the  A/D  limitations  on 
bandwidth.  Ths  sampling  rate  of  the  A/D  needs  only  to  be  twice 
the  frequency  content  of  the  stretched  signal  over  the 
frequency  span  of  interest  ( ie  Af  in  figure  8).  There  is  an 
additional  load  placed  on  the  programable  signal  processor  in 
that  it  is  required  to  perform  a  range  FFT. 

The  actual  design  implementation  of  any  of  the 
previously  stated  spread  spectrum  systems  presents  several 
challenges  and  limitations.  These  problem  areas  generally 
consist  of  hardware  limitations  and  difficulties  associated 
with  spread  spectrum  waveforms. 

State  of  the  art  hardware  limits  the  possible  operation 
of  a  spread  spectrum  system  in  several  areas.  Current  RF 
equipment  restricts  the  bandwidth  to  a  few  hundred  MHz  with 
available  signal  processors  and  A/D  circuitry  limiting  the 
receiver  bandwidth  to  less  than  100  MHz.  Wideband  antenna 
design  becomes  difficult  when  constrained  by  such  parameters 
as  narrow  beamwidth ,  high  gain,  and  low  sidelobes;  all  of 
which  are  important  from  both  an  ECCM  standpoint  and  for 
optimum  target  detection/clutter  suppression  in  a 
fighter/attack  mission.  The  use  of  phased-array  antenna 
systems,  althrough  costly,  provide  the  desired  bandwidth 
response  necessary  for  a  spread  spectrum  radar. 

Another  area  of  concern  is  in  the  proper  use  and/or 
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design  of  wideband  elements.  The  important  figures  of  merit 
include:  phase  and  amplitude  distortion,  RF  power  levels, 
dynamic  range  and  noise  figure.  Most  of  the  phase  and 
amplitude  distortion  is  incountered  in  the  transmission  lines 
which  may  be  coaxial  cable  or  waveguide.  The  frequency 
response  for  long  coaxial  lines  is  generally  non-uniform 
resulting  in  amplitude  distortion.  Properly  matched  lengths  of 
up  to  100  feet  can  achieve  bandwidths  of  approximately  25%. 
Radars  operating  at  or  above  S-band  utilize  waveguide  in  order 
to  avoid  propagation  losses  of  the  higher  modes.  However, 
wideband  operations  create  dispersive  losses  resulting  from 
the  propagation  of  the  TM  modes.  The  delays  are  not  constant 
with  frequency.  These  distortions  can  be  corrected  via 
waveguide  filters,  proper  selection  of  waveguide  size,  and 
digital  correction  in  the  radar  signal  processor.  As  far  as 
power  amplifiers  are  concerned,  the  traveling  wave  tube  can 
achieve  bandwidths  of  10-30%  and  is  generally  the  transmitter 
type  found  in  coherent  airborne  radar.  Also,  the  use  of  solid 
state  devices  such  as  GaAs  FET  amplifiers  provide  low  noise 
(noise  figure  3dB),  high  dynamic  range,  wideband,  and  low 
distortion,  RF  receivers. 

Spread  spectrum  modulation  can  be  classified  as 
intrapulse,  interpulse,  or  hybrid  as  was  previously 
illustrated  in  table  I.  The  technological  limitations  of 
using  these  modulations  is  the  subject  of  this  section.  The 
pseudo-random  polyphase  intrapulse  codes  are  common  in  both 
radar  and  communications  applications.  The  technology  is  well 


developed.  The  hardware  is  readily  available,  low  cost, 
standard,  and  broadband.  Polyphase  code  generators  operating 
at  or  above  500  MHz  are  currently  available. 

On  the  receiving  side,  despreading  processors  or  pulse 
compressors  (PC)  can  be  either  correlation  or  convolutional 
matched  filter  processing.  A  matched  filter  convolves  the 
return  with  the  filters  impulse  response  (a  time  inverted  and 
shifted  form  of  the  transmitted  signal).  The  output  of  the 
filter  is  a  narrow,  time  compressed  pulse.  State  of  the  art 
technology  available  to  perform  the  operation  include: 

1)  digital:  bipolar , CMOS/SOS ,  charged- 
coupled  devices  (CCD) 

2)  discrete  time  analog  CCDS 

3)  surface  acoustic  wave  devices  (SAWD): 
fixed-code  and  programable, coherent 
memory , convolvers ,MOS-hybrid . 

Thus,  current  technology  exceeds  present  or  near  term  system 
requirements,  capabilities,  and  frequency  allocation 
restrictions  for  matched  filters. 

The  generation  of  either  interpulse  or  intrapulse 
frequency  hopped  (FH)  signals  requires  a  fast,  coherent 
switching  frequency  synthesizer.  Two  methods  exist  for 
generating  a  FH  waveform:  direct  and  indirect  synthesis. 
Direct  synthesis  produces  the  output  frequency  from  one 
precise  frequency  source  and  uses  mathematical  multipliers, 
dividers,  adders,  and/or  subtractors.  This  allows  faster 
switching  times  and  fine  frequency  resolution  assuming  the 
time  delay  through  the  bandpass  filters  is  negligible  compared 
to  the  switching  transient.  Indirect  synthesis  yields  the 


output  frequency  by  mixing  one  or  more  voltage  controlled 
oscillators  which  are  often  phase  locked  to  a  frequency 
reference.  In  either  of  the  above  cases,  the  ability  to 
produce  the  frequency  translations  of  FH  systems  is  presently 
available. 

Another  hardware  limitation  involves  the  radar  digital 
signal  processor.  Its  function  is  to  condense  the  sensor  data 
several  orders  of  magnitude  so  that  the  radar  data  processor 
(RDP)  can  make  decisions  concerning  the  radars  current 
tasking.  Figure  9  is  a  generalized  diagram  of  a  typical 
coherent  airborne  radar.  Typical  values  for  a  conventional 
radar  would  include  a  3  MHz  operating  rate  for  the  A/D 
converter  at  the  signal  processors  front  end.  For  a  spread 
spectrum  radar,  the  A/D  may  be  required  to  operate  at  rates  up 
to  the  signal  bandwidth.  This  will  depend  on  the  particular 
system  configuration  as  was  reviewed  previously.  Present 
capabilities  of  A/D  converters  have  been  greatly  enhanced  by 
the  integrated  circuit  chip.  Current  monolythic  A/D  converters 
are  operating  at  several  MHz  with  12  to  14  bits  of  resolution. 
The  tactical  f ighter/attack  mission  requires  this  level  of 
resolution  due  to  the  enviroment  in  which  its  operates  ( ie  ECM 
and  clutter).  Sampling  rates  for  intrapulse  phase  codes 
require  the  A/D  to  operate  on  the  order  of  the  transmit 
bandwidth,  which  could  be  as  much  as  500  MHz.  FH 
configurations,  however,  can  be  designed  to  require  much 
lower,  realizable  sampling  rates.  Future  predictions  for  A/D 
converters  are  promising  in  terms  of  their  effect  on  spread 


spectrum  technology.  GaAs  FET  devices  are  scheduled  to  deliver 
rates  above  1  GHz  with  12  to  14  bits  os  resolution  by  1985 
while  Josephson  Junction  technology  will  provide  10  GHz  by 
1990. 

The  quantized  data  output  of  the  A/D  is  delivered  to 
the  radar  signal  processor  where  algorithms  perform  various 
inner-products  on  large  data  sets.  Such  operations  include: 
filtering,  convolution,  correlation,  and  spectrum  analysis 
(discrete  and  Fast  Fourier  transforms).  Generally,  the 
processors  are  programmable  due  to  the  complexity  ( ie  number 
of  operational  modes)  of  state  of  the  art  radar  systems.  The 
use  of  spread  spectrum  waveforms  places  considerable  burden  on 
the  programmable  signal  processor ( PSP ) ,  where  a  500  MHz  signal 
bandwidth  may  require  the  PSP  to  operate  at  1  giga-operation 
per  second(GOPS) .  This  is  currently  unrealizable  since  present 
PSP  capability  is  in  the  area  of  about  20  MOPS.  Projected 
VLSI,  GaAs  technology  should  yield  1  GOPS  within  10  years.  The 
conclusion  to  draw  from  the  above  discussions  of  the  various 
hardware  limitations,  design  challenges,  and  future  technology 
is  that  as  far  as  device  technology  is  concerned,  spread 
spectrum  is  a  viable  concept  for  coherent  airborne  radar  in 
the  near  future  (10  years).  The  development  and  implementation 
of  spread  spectrum  radar  with  current  technology  is  feasible, 
but  only  if  certain  problems  (such  as  A/D  sampling  rates)  can 
be  circumvented.  Also,  the  transmit  waveform  must  be  carefully 
selected  as  it  has  direct  bearing  on  the  receiver 
configuration  and  thus  feasibility  of  the  design. 


Certain  problems  unique  to  excessively  broadband 
waveforms  exist  that  are  not  encountered  with  low  resolution 
waveforms.  The  following  discussion  presents  the  major  design 
challenges  of  selecting  a  suitable  spread  spectrum  waveform. 
These  characteristics  are  dependent  solely  on  the  waveform 
modulation  (FH, biphase  coded)  structure  and  its  corresponding 
spectral  characteristics. 

Waveform  self-clutter  is  a  problem  with  any  radar 
signal  and  can  be  thought  of  as  any  output  from  a  matched 
filter  due  to  non-matched  targets  ( ie  targats  with  doppler  and 
delay  which  should  not  cause  a  response) .  The  figure  of  merit 
which  is  used  to  describe  the  self-clutter  of  a  waveform  is 
the  ambiguity  function  and  its  corresponding  ambiguity 
diagram.  The  function  is  defined  as: 

"X-CJ-jV)  =  J  s{  t)s*(  t-“V)exp[  j2nvt]  dt  (4) 

The  function  is  the  correlation  of  the  signal  with  a  doppler 
shifted,  range  delayed  version  of  itself.  By  taking  the  square 
of  the  functions  magnitude  and  plotting  it  against  delay(T  ) 
and  doppler(v),  we  obtain  a  visual  representation  of  the 
output  of  a  matched  filter  matched  to  the  signal  and  whose 
input  is  any  combination  of  delayed  and/or  doppler  shifted 
versions  of  the  signal.  This  can  also  be  viewed  as  the  output 
of  a  correlation  followed  by  envelope  detection.  A 
cross-section  of  the  ambiguity  surface  taken  parallel  to  the 
doppler  axis  and  at  T  =  Y0 would  represent  the  filter  output 


corresponding  to  a  fixed  range  (3£)  and  doppler  mismatched  over 
all  possible  values.  The  radar  receiver  approximates  the 
ambiguity  function  via  a  two-dimensional  array  of  correlators 
(time  delay  .vs.  doppler  shift)  or  a  single  array  of  matched 
filters  (doppler  shifted).  In  order  to  resolve  a  target  in 
both  range  and  doppler,  the  majority  of  the  filter  response 
must  be  located  at  the  center  of  the  ambiguity  diagram.  This 
is  where  the  T/v  combination  for  that  target  intersect. 
Therefore,  the  ideal  ambiguity  diagram  would  be  a  single  spike 
at  the  origin  and  zero  elsewhere.  This  is  physically 
impossible.  The  reason  for  this  is  made  clear  in  Chapter  III, 
where  a  mathematical  treatment  of  the  ambiguity  surface  is 
presented. 

A  receiver  response  to  interference  ( clutter , noise , non- 
matched  targets)  is  largely  a  function  of  its  ambiguity 
characteristics.  Complicating  the  problem  is  the  fact  that 
there  does  not  appear  to  be  an  analytic  method  to  construct  a 
signal  from  a  particular  ambiguity  surface.  Therefore,  a 
particular  signal (s)  is  selected  and  from  this  we  can 
construct  the  ambiguity  surface  and  predict  the  waveforms 
self-clutter  properties. 

The  inherently  high  range  resolution  of  a  spread 
spectrum  radar  (assuming  the  entire  bandwidth  is  effectively 
processed)  presents  a  problem  termed  "target  super¬ 
resolution".  This  super  resolution  is  a  direct  consequence  of 
the  Fourier  integral  property:  large  spectral  extent 
transforms  to  short  time  extent.  For  a  target  whose  spatia- 


extent  is  greater  than  the  radar  resolution  (waveform 
bandwidth),  it  appears  to  the  receiver  as  several  targets  each 
resulting  from  the  targets  individual  scatterers.  Assumming  a 
transmit  bandwidth  of  several  MHz,  the  resolution  is 
approximately  one  foot.  Thus,  a  forty  foot  target  is  "broken 
up"  into  forty  seperate  returns.  If  the  radar  receiver 
bandwidth  were  matched  to  the  largest  target  extent,  the 
individual  returns  could  simply  be  non-coherently  integrated 
and  the  complete  RCS  recovered.  The  individual  scatterers  of 
the  super-resolved  target  may  be  too  weak  to  be  detected.  This 
is  the  essence  of  the  super-resolution  problem  and  the  focal 
point  of  this  paper.  In  order  to  compensate  for  the  over 
resolution  properties  of  wideband  radar  signals,  a  "target 
filter"  is  required.  This  filter  will  perform  an  integration 
on  the  individual  returns  in  one  of  two  ways:  post-detection, 
non-coherent  range  collapse  or  pre-detection,  coherent 
integration. 

The  non-coherent  target  filter  consists  of  post¬ 
detection  integration  or  equivalently  range  collapsing  the 
range  data  to  match  the  target  dimension.  This  is  done  after 
envelope  detection  and  requires  the  entire  signal  bandwidth  to 
be  processed.  If  the  majority  of  the  target  scatterers 
occupied  a  single  range  cell, then  the  optimum  processing  would 
be  single  cell  observation.  Target  scatterers  which  have  low 
RCS  and  occur  in  several  range  cells  require  cell-to-cell 
integration  for  detection  enhancement.  Processing  the  entire 
signal  bandwidth  places  a  burden  on  current  programmable 
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signal  processor  capabilities.  In  certain  instances,  500  MHz 
bandwidths  with  phase  coding  for  example,  this  makes  the 
system  unrealizable.  Other  disadvantages  of  post-detection 
integration  include  those  normally  associated  with  this  type 
of  processing:  collapsing  loss  and  post-detection  integration 
loss.  In  an  effort  to  reduce  the  bandwidth  restrictiions  and 
PSP  loading,  a  coherent,  pre-detection  filter  is  considered. 

A  coherent  target  filter  is  a  matched  filter  ( ie 
correlator)  which  may  depend  on  the  return  signal  being 
narrowband.  Assuming  the  signal  spectrum  is  represented  as 
S(jw),  and  the  spectral  characteristics  of  the  scattering 
function  (range-extended  target)  as  H(jw),  a  matched  filter 
would  have  a  spectrum  which  is  S ( - jw )H ( - j w ) .  The  bandwidth  of 
this  filter  may  be  narrowband  (see  Chapter  III)  with  respect 
to  the  signal  spectrum.  If  this  conjecture  holds,  a 
considerable  processing  advantage  (on  the  PSP)  could  be 
gained;  the  result  being  realization  of  a  spread  spectrum 
radar  system. 

The  objective  of  this  thesis  will  be  to  explore 
narrowbanding  effects  of  the  transmitted  waveform.  This  will 
be  accomplished  via  comparisons  of  signal -to-noi se  ratios  for 
non-coherent  filtering  .vs.  coherent  narrowband  filtering. 
Two  waveforms  are  considered:  FH  and  biphase  coded.  Also, 
target  models  will  include  Swerling  I  and  II  amplitude 
fluctions.  The  thesis  will,  in  the  process  of  simulating  the 
above  filter,  attempt  to  verify  the  results  of  Salzman  (Ref'll) 
that  the  coherent  filter  is  advantageous  (SNR)  only  20%  of  the 


Ill  Theoretical  Considerations 


Complex  Signal  Representation 

Many  radar  signals,  systems,  and  filters  of  practical 
interest  satisfy  the  narrowband  criterion.  This  means  the 
frequency  content  (signal)  or  frequency  response  (filters)  is 
confined  to  a  frequency  band  which  is  a  small  percentage  of 
the  center  or  carrier  frequency  (Ref 12: 55).  A  generally 
accepted  figure  of  10%  has  been  applied  in  most  radar 
applications  and  will  be  adhered  to  through  out  this  study 
(Ref2).  Complex  signal  notation  will  now  be  presented  assuming 
narrowband  signals.  This  is  not  a  necessary  condition  but 
serves  to  simplify  interpretation  of  the  results.  A  real, 
narrowband  signal,  f(t),  can  be  expressed  as 

f(t)=  Re[u( t)exp( jwG  t) ]  (5) 

Since  u(t)  is  by  definition  a  slowly  varying  function  of  time, 
its  Fourier  transform,  U(w),  is  band  limited  to  w^-wQ .  The 
exponential  factor  serves  to  shift  the  spectrum  to  the  right 
by  w0  ,  resulting  in  a  purely  positive  spectrum.  The 
component,  u(t),  is  referred  to  as  the  complex  envelope  of  the 
real  signal,  f(t),  and  contains  the  angle  and  amplitude 
modulation  of  the  real  signal  (Ref  3: 58).  The  product 
u(t)exp(jw0  t)  represents  the  analytic  signal  and  is  further 
developed: 

f(t)=  Re[u( t)exp( jw  t) 1  (6) 


(8) 


where  fB(t)=  f(t)+jf(t) 

f(t)  is  the  Hilbert  transform  of  f(t)  and  is  defined  by 

+  0» 

f (t)=  i/n  r  s<fr>  d?* 

where  the  integration  calculates  the  Cauchy  principle  value. 
The  analytic  signal  spectrum  has  only  postive  frequencies  and 
is  related  to  the  real  signal  spectrum  via 

U(w)  =  (2F(w)  ,  w>0 
J  F ( W )  ,  w=0 
(  0  ,  w<0 

u(t)  is  a  complex,  lowpass  function  which  can  further  be 
written  as 


u(t)=  x(t)+j  y ( t ) 

where  x(t)  and  y(t)  are  the  quadrature  components  of  f(t).  In 
this  manner  we  can  mathematically  describe  signals  in  the 
complex  domain  and  derive  the  physical  process  via  the  above 
transforms.  Further,  any  bandpass  signal  can  be  mathematically 
processed  in  complex  signal  notation  and  physically 
interpreted  as  manipulations  of  the  quadrature  components 
(Ref2).  To  illustrate  this  we  start  with  equation  7  and  a 
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bandpass,  real  signal  Z(t): 


Z(t)=  Re[Za ( t) ] 

=  Re[u(t)exp( jw0 t) ] 

=  Re[  ju(t)|  exp(  jwQ  t)  exp(j$t)] 


where  <t>(t)  is  the  angle  of  u(t). 


Z(t)=  Ret^X2  <  t)+Y2  ( t)  exp[w0  t+<J>(  t )  ] 
=  ^X*(  t )  +Y*(  t )  cos[w0t+ <t>(  t)  ] 


(9) 


using  Z(t)=  X(t)cos(w0t)-Y(t)sin(w0t) 

=  M(t)cos{wot+0(t) )  (10) 


where  M(  t )  =  V  ( t )  +Ya  ( t ) 

0(t)=  tan1fY(t)1 

IxTtlJ 


and  comparing  this  to  equation  9  we  see  that  Z(t)  can  be 
represented  as 


Z  ( t )  =  Re(  u(t)  exp(jwat)  exp(  j<j>(  t ) )  ] 

=  Re[a(t)  expj<t>(t)  expjw0t] 

where  a(t)=  M(t)=  |u(t)| 

and  0(t)=<£(t)  above. 

Thus,  u(t)  contains  all  the  angle  and  amplitude  information  of 
the  real  bandpass  signal.  Suppression  of  the  exp(jw0t)  factor 
causes  no  loss  of  information;  therefore  the  complex  envelope 
u(t)  will  be  used  to  represent  the  signals  simulated  in  this 
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study.  The  above  has  assumed  the  narrowband  criterion  so  that 
a(t)  (the  signal  amplitude  modulation)  can  be  directly 
associated  with  |u(t)|  and  <J>(  t )  associated  with  any  angle 
modulation  (0(t)). 

In  an  analogous  manner,  complex  notation  can  be 
extended  to  the  definition  of  filter  transfer  functions.  The 
real  impulse  response  of  a  linear,  time  invariant  filter  may 
be  expressed  as 


h(t)=  Re[2  h(t)  exptjijt)] 


(11  ) 


where  h(t)  is  the  complex  impulse  response  function  and  the 
factor  of  2  is  for  convenience  in  later  expressions.  Let  x(t) 
and  y(t)  be  the  real  input  and  output  expressions  respectively 
and  express  their  analytic  counterparts  as 


Further, 


xa(t)=  x(t)+  jx(t)=  x(t)expjfit 
ya ( t ) =  y(t)+  jy(t)=  y(t)expjfit 

y( t ) =  h(t)*s(t> 


where  *  denotes  convolution  and 


(12) 

(13) 

(14) 


ya(t)=  h( t ) *x( t ) + jH[ h( t )*x( t ) ) 

=  h(t)*x(t)+jh(t)*£(t)  (15) 


where  H  represents  the  Hilbert  transform,  and  we  have  used  the 
property  of  Hilbert  transform  theory:  if  y ( t ) =v( t ) *x ( t ) ,  then 


(16) 


y ( t ) =v( t ) *x( t ) .  Proceeding  further, 


y  (t)=  h(t)*[x(t)+jx(t)] 
=  h(t)*x.(t) 


Equation  1 6  can  be  written  in  integral  form  as 


also. 


ya  ( t )  =  f h(T)x  (t-T)  dT 
y(t)expjfit=y'h('J')x(t-T)expj  (t-T)  dT 
-  y(t)=yli(T)exp-jfiTx(t-T)  dT 


using  the  property  for  any  complex  function,  z(t), 

Re[z(t)]s  z(t)/2  +  z*(t)/2  and  applying  this  to  2h(T)expj/lt 
we  obtain 


h(T)=  h(T)expjAT  +  h  (T)exp-jiVT 


(18) 


h(T)exp-jAT=  h(T)  +  *  (T) exp- j2hT- 


Equation  18  is  now  substituted  into  equation  17  to  obtain 


y  ( t ) 


= j h(T)x( t-T )  dT  h*(T)x(t-T)exp-j21iTdT 


For  narrowband  signals,  the  second  integral  becomes  negligible 
( ie  designates  radian  frequency)  and 


y(  t)  =  /  h(  T)x(  t-T)  d^ 


We  see  that  the  complex  envelope  of  the  output  is  the 


convolution  of  the  complex  envelope  of  the  input  with  the 
complex  impulse  response  of  the  filter  (Ref12:74).  The  real 
bandpass  output  is  obtained  via 


y(t)=  Ret y( t )expjw0 t ] 

=  Ret (x( t)*n( t) )expjw0 t]  (21) 

Given  the  transfer  function  of  a  narrowband  filter,  the 
complex  impulse  response  is  determined  by  the  following: 

H(jw)=  H(w-A)  +  ( -w-a)  (note  equation  18  above)  (22) 

For  a  filter  with  lumped  circuit  elements,  the  poles  are 
located  near  ft  .  If  the  pulse  near  +  fl.  are  associated  with 
H(w-A)  and  the  poles  near  -  £l  associated  with  H(w+Jl),  then 


and 


h(t)=  Re [2  h(t)expjj(lt] 


H(jw)expjwt  dw 


Therefore,  the  procedure  is  to  take  the  positive  frequency 
content  of  the  real  transfer  function,  shift  it  to  the  origin, 
and  inverse  Fourier  transform  to  determine  h(t)  (Ref  12:74). 
This  can  be  seen  graphically  through  the  following  three 
spectral  representations: 


C.  wwu union  rw^cnow 


^ _ 1 - i - - - 

Figure  10.  Graphical  Representation  of  Complex  Signal  Spectrum 

The  above  procedure  will  be  used  to  develope  complex 
representations  for  several  radar  signals  and  the  simulated 
target  filter. 

Matched  Filter  Theory 

A  standard  engineering  parameter  used  to  measure  the 
effectiveness  of  a  radar  receiver  is  the  signal-to-noise 
ratio.  The  ratio  is  defined  as  the  peak  signal  power  divided 
by  the  average  noise  power.  The  goal  of  any  radar  system  is  to 
maximize  this  ratio  thus  improving  the  ability  to  detect  and 
track  targets.  Signal-to-noise  also  has  an  indirect  affect  on 
the  range  and  doppler  estimation.  The  maximization  of  the 
signal-to-noise  ratio  leads  to  the  definition  of  the  matched 


filter.  Identical  results  can  be  achieved  via  decision  theory, 
specifically  the  optimum  receiver  in  white  Gaussian  noise 
(Ref3:143, 145-149) . 


Assume  an  input  waveform  z(t)=s,  (t)+n,  (t)  where:  s,  (t)  = 
input  signal,  n,(t)=  noise  component;  and  an  output  waveform 
y ( t ) =sc ( t ) +nc ( t ) .  The  objective  is  to  determine  the  linear 
filter,  h(t),  which  maximizes  the  ouput  signal-to-noise  ratio 
at  some  particular  time  t=tm.  The  signal-to-noise  ratio  is 
defined  as: 


sj  (tm) 

S/N=  (23) 

no<tm) 


For  white 
function  is 


Gaussian  noise, 
impulsive: 


we  know  that  the  autocorrelation 


Rn  ('$'>= 


(24) 


£ 


where  N0/2  =  amplitude  of  the  noise  spectral  density,  i(r)  = 
Dirac  delta  function.  Thus, 


I  1 2 

fs,  mh(tm-t)  dr 

s/N  Jtr- — 5 - - 

f n0/2  h  (tm-r>  dr 


(25) 


performing  the  Fourier  transform  of  the  above  and  applying  the 
Schwarz  inequality  leads  to 


H{ w)S.  (w)  exp(jwtm)  dw 


2TT  N0/2  w)j  dw 


2TT  Ne/2 


If  H(w)  is  selected  such  that  the  equality  is  satisified,  then 


H(w)=  2k_S*(w)  exp(-jwtm) 
No 


Therefore,  the  optimum  filter  in  white  Gaussian  noise  is 


h(t)=  Zk,  S  (t  -t) 


Next,  we  repeat  the  derivation  using  complex  signal 
notation.  Again  we  seek  the  filter,  h(t),  which  maximizes  the 
signal- to-noise  ratio  at  t  =  tm  •  The  input  is  x(t)= 
Re[x( t )exp( jw0t ) ] ,  the  filter  h(t)=  Re[ h( t )exp( jw0 t ) ] ,  and  the 
output  signal  is  y(t)=  Re[y( t )exp( jw0 t) ] .  White  Gaussian  noise 
is  assumed  to  pass  through  the  bandpass  filter,  thus  R„ )= 
Re[R~  (T)exp(  jw0'?’)  and  R»('?)=N0/2  6(?").  Equation  28  is  utilized 
to  assert  that 


y  ( t )  =  1  /2  j x#(tm-t+?")x(r)  d^ 

=  R2(t-tm) 

using  the  same  definition  of  signal-to-noise  ratio, 


X(f  )H(f  )exp(2TTft  )  df 

m 


The  Schwarz  inequality  produces 


V  | H  (  f  )J  df  f  |x(f)| 

n0  y|H(f)fdf 


selecting  the  equality  yields  H(  f  )  =X(  f  )exp( -jwt„, )  and  via  the 


Fourier  transform: 


h(t)=  x#(tm-t) 


(32) 


In  the  next  section,  the  optimum  receiver  in  white 
Gaussian  noise  is  derived.  The  result  is  a  matched  filter  and 
we  see  that  Decision  theory  and  the  maximization  of  the 
signal -to-noise  ratio  lead  to  the  same  conclusion. 


Optimum  Receiver  in  White  Gaussian  Noise 

The  determination  of  the  optimum  receiver  will  assume 
that  we  recieve  a  signal  with  known  parameters,  and  apply 
Bayes  detection  criteria  in  the  form  of  a  likelihood  ratio. 
Knowledge  of  the  form  of  the  likelihood  ratio  yields  a  test 
statistic  which  is  used  to  interpret  the  structure  of  the 
optimum  Bayes  receiver (Ref 8 : 1 8 ) . 

First,  sampling  theory  is  applied  to  a  bandlimiteu 


*  vr  4.*  «_**  -  •  *  •  - 


•*',•**  *  •  •' 


ws* 


signal,  s(t),  limited  to  (-wc,wc  ).  The  signal  is  sampled  in 
time  every  1/2wc  seconds  and  is  written  as: 


♦  • 


s ( t )  =  ^  s(k/2fc  ) 


sin[2TTfc  (t-k/2fc  )  ] 
2TTfc  (t-k/2fc  ) 


(33) 


where  fc  =wc  /21T  (Ref  3 :  51 -53  ) . 

It  can  be  shown  that 

T 

r  .  *vr  2 

E=  Jsz(t)  dt  «  1/2fc  £  f4  (k/2£;  )  (34) 

<»  K»t 

where  T  is  the  signal  duration.  Equation  34  simply  states  that 
by  sampling  at  twice  the  bandwidth  (Nyquist  rate);  we  can 
represent  the  majority  of  the  energy  in  the  sampled  domain. 

Assuming  a  known  signal  bandlimited  to  (-fc,fc  ),  the 
noise  is  also  bandlimited  and  the  input  waveform  z(t)  becomes 
z,  where  2=5+13  (zj=S(  +n(  )  and  each  sample  z  =z(  i / 2 fc  ), 
s  =s(  i  /2fc  ) ,  and  n  =(i/2fc),  =  1,2,3,....2fcT.  The  hypothesis 

of  the  signal  plus  noise  (Ht)  is  described  by  the  conditional 
joint  probability  density  function 


p(s/Hi  ) 


nm  N  exp(-1/2fcN0  V(z,  -s.  f) 

(  2TT)*  ( N0  f c  V 


(  35  ) 


Where  N  =  2fcT  and  the  power  spectral  density  of  the  noise  is 
No /2 .  In  a  similiar  fashion,  the  noise  alone  hypothesis  ( Ho ) 
is  represented  via 


p(  z/Ho  )  = 


(  2TT)H*(  fe  No  )N 


exp(-1/2ftN.^z*  )  (36) 


The  likelihood  ratio  is 


A(z)  = 


p(z/H1)  exp[ -1/2f*Nc  £  (zi-sL  ft 

"  *  '  ■  ■ 

P(z/H0)  exp(-1/2ftN0  £  z*  ) 


Applying  equation  34  to  the  likelihood  ratio  we  have 

T  7 

A(z)=l(z(t)  )  =exp[  -1  /N*  Jilt)  dt  +  2/N„ /z(  t)s(  t)  dt] 


Further, 


In  l(z(t))=  -E/No  +  2/No  Jz(  t)s(  t)  dt 


Thus,  the  Bayes  likelihood  ratio  receiver  is 


r 

2/N, /z ( t ) s( t )  dt  <  Threshold  +  E/N0 
J  H 


The  left  side  of  the  inequality  or  test  statistic  is  the 
cross-correlation  of  z(t)  and  s(t).  A  matched  filter  is  used 
to  perform  the  correlation  and  we  note  that  the  optimum 
detector  in  white  Gaussian  noise  is  the  matched  filter. 
Specifically,  for  a  fixed  probability  of  false  alarm,  the 
matched  filter  is  the  receiver  that  maximizes  the  probability 


of  detection  (Ref2). 


As  a  final  note  on  matched  filters,  we  observe  that 
given  a  complex  signal  (analytic)  ^(t),  the  response  of  the 
matched  filter  is 


p(V=j  i(t)  i(t-T)  dt 

The  above  complex  autocorrelation  function  can  be  equivalently 
represented  as  a  convolution  (designated  by 

*(**)  *  **(•*)  (4i  > 

Equation  41  will  later  be  used  to  simulate  the  complex  matched 
filter  processing  as  convolutions  ( ie  FFT  algorithm). 

Received  Signal  Representation 

The  basic  radar  signal  examined  in  this  study  is  a 
high-resolution,  coherent  pulse  train  with  some  form  of 
frequency  or  phase  modulation.  The  simulation  will  assume  we 
are  processing  the  video  (baseband)  signal  which  resu  *  frri 
an  enviroment  (target  model)  of  many  point  scatterers 
Simulating  the  received  signal  from  a  single  scatterer  will 
permit  the  construction  of  the  return  from  a  conglomerate  of 
scatterers  via  superposition.  This  will  automatically  simulate 
the  phasor  interference  among  several  scatterers  and  is  the 
only  way  to  effectively  simulate  the  enviroment  of  hiqh- 
resolution  radar  (Ref4:10). 


The  general  form  of  the  received  signal  from  a  point 
scatterer,  moving,  fluctuating  in  amplitude  or  phase,  and 


antenna  scanned  is 


sR(t)=  sT(t-*(t))(  A*7(4n)3  f*)G(t)7(t)  (42) 

where  G(t)  is  the  one-way  antenna  power  gain,  7(t)  is  the 
complex  reflection  coef f icient ( a  voltage  quantity  related  to 
radar  cross-section(  a  ^vf),  and  l*(t)  is  the  round  trip  delay 
for  the  target  (Ref4:6-10).  Equation  42  is  difficult  to 
simulate  and  therefore  several  simplifications  are  necessary. 
If  G(t),  ^(t),  and  7 (t)  are  slow  varying  functions  of  time, 
equation  42  becomes  a  delayed  and  doppler-shif ted  replica  of 
the  transmitted  signal: 


sR(t)=  st(  t-T)exp(  j2TTvt )  (  £  /  (  43TT)3  r4  )G  Y  (43) 


(44) 

(45) 

with  r=  range 

v  =  doppler  frequency  shift 
r=  velocity 
X=  wavelength 


where  ?=2r/c 
v=-2r/\ 


The  above  specifications  are  generally  termed  constant-doppler 
theory  (slowly  moving  target  with  respect  to  a  radar  look 
interval)  and  is  valid  if  the  following  hold: 

1 .  T is  constant  for  r  < <  c/2TB  (46) 


where  TB  is  th  time -bandwidth  product 
and  c=3X10  m/s 


2.  v  is  constant  for  r  << 


/2T* 


(47) 


A  factor  of  four  is  usually  sufficient  to  ensure  the 
conditions  are  satisified  (Ref4:13).  Continuing  with  equation 
42  and  setting 


V.  =  [  X  /  (  4  [  1  /r*  ]G  Y  (  48  ) 

K  k  k  k 

where  the  subscript  k  refers  to  the  k*h  scatterer ,  we  now  have  a 
received  signal  of  the  form 


sR(t)  =  VksT  (t-T(t)) 


(49) 


where  V ( t )  is  not  approximated  as  constant.  In  terms  of 
complex  notation  we  can  express 


sT(t)=  uT( t)exp( jwc t)  (50) 

and 

s„<t)=  V„uT[t-m)]exp[j2Hfe(t-*<t))]  (51) 

The  delay  is  a  function  of  time  and  depends  on  the  target 
range  at  the  instant  of  signal  reflection: 


>(t)=  2/c  r[t-^t)/2] 


(  52) 


-  -  ."V  -s.  ’ 


V’l  *«>  V  *.  V.  .  •*_  *. 


where  r(t)  is  the  range  function.  Performing  a  Taylor  series 
yields 


*?(  t )  =  'V,,  -vu  ( t-tp  )  +  additional  terms  (53) 

yK  ,and  vK  represent  the  time  and  doppler  coefficients  of  the 
scatterer,  and  to  is  the  time  at  which  they  are  defined.  If 
equations  46  and  47  are  satisified  we  can  write 

sR(t)=  V*  uT  ( t  -  %  )  expj  2TT[  fc  ( t- ^k.)  +vK  ( t-t„  )  ]  (54) 

Further,  we  assume  multiple  scatterers,  thus 

sR(t)=  exp(  j2Hfc  t)^^VRuT(t-Tt)exp(  j2FvKt)  (55) 

K 

The  above  equation  will  be  simulated  with  the  carrier  term 
suppressed  and  a  simplification  made  which  applies  to  noise- 
type  waveforms. 

If  the  amplitude  or  phase  modulation  of  a  radar  signal 
is  noise-like,  equation  55  cannot  be  simplified  further  unless 
we  assume  short  pulses.  For  a  pulse  train  as  shown  in  figure 
12,  the  bandwidth  is  approximately  the  inverse  of  the  pulse 
width  B=  1/Tp  ,  and  the  time-bandwidth  product  is  on  the  order 
of  unity  (TB=1 ) . 

Resolution  in  delay  is  given  bydT  =  1/B  and  is  also 
determined  by  the  pulse  width:  ^?=Tp.  When  Tp  is  short  enough 
to  achieve  a  given  resolution  in  delay,  it  will  be  impossible 


46 


Figure  11.  Coherent  Pulse  Train  (Ref 4: 55) 

to  resolve  scatterers  in  doppler  due  to  the  property  Av=1/Tp  , 
thus  it  is  possible  to  assume  all  scatterers  have  the  same 
doppler  velocity  (Ref4:51).  Further.,  the  doppler  is  assumed 
equal  to  zero  for  convenience.  This  yields  a  received  signal: 

sR(t)=  exp[  j2TTfc  t  ]  (56) 

K 

Equation  56  (with  exp[j2¥fct]  suppressed)  will  be  directly 
simulated  in  this  study.  The  assumption  of  short  pulses  (TB=1 ) 
is  equivalent  to  assuming  the  doppler  spread  of  all  scatterers 
is  small  (Ref4:52).  It  is  the  later  assumption  which  we  are 
simulating.  The  actual  waveform  may  or  may  not  have  a  time- 
bandwidth  product  on  the  order  of  unity. 

Ambiguity  Function 

In  order  to  quantitatively  assess  the  performance  of  a 
particular  radar  waveform  and  its  procersinq  filter(s),  we 


need  a  function  which  describes  the  signals  performance  in  a 


multitarget  enviroment.  This  function  is  the  ambiguity 
function  and  is  a  measure  of  the  signals  range  and  doppler 
resolution  (Refill  27).  Resolution  is  the  ability  of  the  radar 
to  resolve  (detect)  two  or  more  targets.  For  the  purposes  of 
this  study,  the  multitarget  enviroment  consists  of  point 
targets  located  on  the  same  target.  Further,  The  ambiguity 
function  can  be  though  to  represent  the  self-clutter  of  the 
waveform  as  well  a  measure  of  performance  in  a  multitarget 
scenario.  This  will  become  evident  as  the  properties  of  the 
function  are  stated. 

Suppose  we  have  two  point  scatterers  with  different 
ranges  and  we  are  using  matched  filter  processing.  Let  x(t).= 
a(t)cos[w0  t+<Mt)J  be  the  return  from  scattererl  and  x(t-T)  = 
a(  t-Y)cos[w;(  t-?)+  <t>  (t-T)  ]  be  the  return  from  scatterer2.  The 
range  seperation  between  the  two  is  cl'/2,  where  c  is  the  speed 
of  light.  We  will  next  use  complex  signal  representations: 
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x(t)=  Re(u(t)exp( jw« t) ] 


or 


?=  y*Jx(t)-x(t->)|  dt 

e= J [u(  t)-u(  t-Y)exp( -jw0?*)  ]  [u#(  t ) -u*(  t- f)exp(  jw.»  ]  dt 


(57) 


48 


=2J  |.(t|  -2  Re[exp( -jw0T")  u(  t  )u(  t-'J')]  dt 

="*2[  2E]  -  2  Re[B(»exp(-jw0l-)  ]  (58) 


B(  *}■  )  = 


/ 


u(t:u(  t-r) 


maximum  amount  of 
parameter : 


dt.  Minimizing  B(h)  will  result  in  the 
resolution.  Defining  a  delay  resolution 


[\b^)\  dr 

/jM(f)|2df 

|B(0)f 

(  /|M(f)|a  df)‘ 

(59) 


Equation  59  utilizes  Fourier  transform  properties.  Further,  a 
range  resolution  constant  can  be  defined  as  R=cTr/2  and  an 
effective  bandwidth  constant  as  B#  =1 /2TR  .  The  point  is  that 
improving  range  resolution  requires  increased  effective 
bandwidth. 

By  analogy,  we  can  determine  the  velocity  resolution  of 
two  point  targets  where  x(t)  is  the  first  target  return  and 
x(t)exp(jwD  t)  is  the  doppler-shif ted  return  of  the  second 

— A 

target.  Again  e  is  formed  and  we  proceed  as  before: 

i 


= y~|x(  t)-x(t)exp(  jwQ  t)|  dt 

=  2^*|x(  t)|*  dt  ~Re^J 


(60) 


-Re [ I  x(t)x(t)exp( jwQt)  dt]  (61) 


using  the  transform 

t)x*  ( t  )exp(  j21Tvt )  dt  «— *^X*(v)X(v+w 
=  2 TTV  and  X(v)=Fourier  transform  of  x(t).  Next,  we  define  the 


)  dw  =K.(v)  with  w_ 


frequency  resolution  constant  T,  as 


(62) 


which  leads  to  the  effective  time  duration  =  1/Ff  .  Now  for 
two  scatterers  at  the  same  range,  they  can  be  resolved  in 
doppler  if 


A  wd/2  >  Fr 


where  is  the  doppler-shif t  difference. 

A  similiar  consideration  of  target  resolution  for  a 
single  parameter  (range  or  doppler)  is  related  to  the  physical 
processing.  Return  signals  are  effectively  filtered  in  bins  of 
range  and  doppler.  A  bin  width  or  interval  is  defined  as  the 
half -power  width  of  the  central  lobe  of  the  matched  filter 
response.  Restating  B(T)=J' u*  ( t  )u(  t-t)  dt  =J^ |m(  f  )jexp(  j 2TTf t )  df 
and  defining  the  bin  width  as 


assuming  B(0)=1,  then  1  /2  =  1  +1  /2( -2$')?" 
where  f  is  the  rms  bandwidth  and 
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a 


'  w  2  F ( w )  dw 
/  |F(w)|*  dv 


This  yields  =1/2  andV  =  1//2£  .  The  response  envelope 

spectrum  is  two-sided,  therefore  we  use  T  =  2V  =  ■17/$  and  define 
the  range  bin  as  rfe  =Vc/2=c/ jT^.  For  velocity  resolution,  we 
note  that  R^v)]  =jfju(t)[Z  exp(jznvt)  dtj  .  Again,  performing  a 
Taylor  series  expansion  about  v=0,  we  obtain 

K  (V)  =  1/2  =  Ik  (Oil  .I/^/Ik  ( v >[  .  V*  (66) 

“  ~ j;4 

Assuming  normalized  signal  energy,  Ko(0)=1,  we  define  V^sJT/tj 
where  tj  is  the  rms  time  duration  defined  as 

2  /t8k(t)|ldt 

td=77 — (67) 

/  jf  ( t  )lcdt 


The  doppler  bin  width  is  r^  =  AV^/2  =>//2’tj  (68) 

The  above  can  be  extended  to  the  more  general  case  of 
two  scatterers  which  differ  in  both  range  and  doppler.  In 
similiar  fashion  the  integral  square  error  is  maximized; 


/|x(t)exp(jw,t)-x(t-*')l  dt 

r  i*  "r 

=**2  /|x(  t )  dt  -  Re[j  x(t)x{t-T)exp(  jwet)  dt]  (69) 


Maximizing  e  requires  the  minimization  of 


x(  t)x(  t-»exp(  jwDt)  dt 


Appling  complex  envelope  notation  results  in 


u(t)u*(t-‘J')exp( jw0t)  dt 


(70) 


Equation  70  is  the  ambiguity  function  of  the  waveform.  The 
ambiguity  function  is  the  complex  envelope  of  the  matched 
filter  response  in  delay  and  doppler  (Ref9:112).  The  real 
response  is  calculated  by  taking  the  real  part  of  X(  X d  ) 
exp(jwot)  (Ref9:112).  Ifj^'y(T,w)|  is  calculated  and  plotted  in 
the  X plane,  a  3-dimensional  diagram  is  generated.  The 
ambiguity  surface  [  j‘X(?T w )J  ]  visually  depicts  responses  of  the 
matched  filter  for  various  ?,w  combinations.  The  central  peak 
of  the  diagram  (?’=0,v=0)  is  the  filter  response  for  a  perfect 
match,  and  all  other  response  peaks  are  resultant  from 
scatterers  (or  targets)  which  are  non-matched.  Thus,  any 
waveform  which  generates  sufficiently  high  subsidiary 
responses  has  potential  targets  at  the  corresponding  -$\w 
values.  In  the  case  of  a  single  target,  more  than  one  target 
might  be  incorrectly  detected  or  in  the  case  of  several 
targets,  only  one  might  appear  to  dominate.  The  ideal 
ambiguity  surface  would  be  an  impulse  located  at  T=v=0 
( Ref 9 :113— 114) .  The  objective  would  logically  be  to  generate  a 
waveform  which  has  this  characteristic.  Unfortunately, 


specification  of  the  ambiguity  function  does  not  uniquely 
define  a  waveform.  The  utility  of  the  ambiguity  diagram  can  be 


further  understood  by  developing  several  properties  of  the 
function. 

In  the  frequency  domain  we  have: 


y  (T/W)=y*M(f+w)M(f  )exp(-j2TTft)  df 


(71  ) 


and  we  note  that  a  maximum  occurs  at  the  orgin 


'fcftw)  <  Y(0,0)  =  2E 


Also, 


=  *(T,w) 


"Yl  "^w)!  =  |Y<r,w)|  therefore  Y(-T,-w) 
function  is  symetric  about  the  orgin  (Ref9:1 18-126) . 

Next,  we  explore  the  volume  under  the  function: 


and  the 


volume= 


d^dw 


If  we  take  the  two-dimensional  Fourier  transform, 


■  *,v  t  |V<  r,  v  »|)  =^Y<  r,  v  ilexp  ( - 


j  )exp(  -  j  2TTuv  )  d?dv  = 


(72) 


We  see  that  the  ambiguity  surface  is  its  own  two-dimensional 
transform  and  further  by  letting  u=^=0  we  have: 


"X(^v) 


2 


d?"dv 


(73) 


Equation  73  states  that  the  area  under  the  surface  (and  thus 


s  *. 


the  ambiguous  responses)  is  constant  and  equal  to  4  (assuming 
normalized  signal  energy, E=1 )  Thus,  ambiguous  returns  may  be 
shifted  within  the  \ ,v  plane  but  the  sum  total  remains 
constant  (Ref 9: 1 20-1 22 ) .  The  objective  is  to  design  a  waveform 
whose  ambiguity  diagram  has  the  response  peaks  distributed  to 
minimize  the  self-clutter .  Targets  at  ranges  and  velocities 


such 


that 


Y<  t  /V) 


IS 


V  ( o ,  o )  | 


are 


comparable  to 

indistinguishable  to  the  radar  (Ref8:26).  The  width  of  the 
response  peak  aboutj  ^(  0 , 0  )|  defines  the  waveform  resolution 
while  the  peaks  located  away  from  the  orgin  are  ambiguous 
(Ref8 : 26 ) . 


Sampled  Signal  Theory 

In  any  digital  computer  simulation  the  radar  signals 
must  be  sampled.  This  section  derives  the  Fourier  transform 
relationships  of  sampling  in  the  time  domain.  Frequency  domain 
sampling  is  derived  in  analogous  manner. 

Consider  an  arbitrary  complex  signal,  u(t),  with 
Fourier  transform  U(f)  and  energy  Ec  .  Sampling  u(t)  at 
discrete  times  kAt  yields  an  approximation  of  U(f): 

Us(f)=  At  £u(  kAt )  exp(  -  j2itf  kat )  (74) 

K'-«e 

A  function  sampled  in  the  time  domain  will  be  repetitive  in 
the  frequency  domain  (Ref 4:85).  Thus, 

-♦oo 

Us(f+m/At)=  At  2^u(  kAt )  exp- j  2  rr(  f  k  At +mk  )  (75) 

K*.co 
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The  product  mk  is  an  integer  which  leaves  the  phase  unchanged 
and  the  mk  term  is  omitted  (Ref9:85).  U§(f)  is  a  periodic 
function  with  frequency  ff  =1/  At.  Replacing  the  Fourier 
transform  for  u(kAt)  in  equation  74  we  find: 

4e» 

(76) 


r 

Us(f)=  I  U(  f  ^ )  [  At^^exp- j  2"n  ( f -f  * )  kat  ]df / 

— '  Q0  ||*S  gQ 

and  noting  that 


we  obtain 


K 


At/  pxp-j2TtfAt  =  (  f-mfr  ) 

4  ® 

Us(f  )  =  2U<f-mfr  ) 


(77) 


-  <30 


which  states  that  the  spectrum  due  to  time  sampling  is  the  sum 
of  the  original  continous  spectrum  shifted  at  intervals  equal 


Figure  12.  Spectra  of  Continous  and  Sampled  Signals  (Ref4:86) 


The  two  spectra  are  comparable  for  fr=1/At,  and  alaising  is 
minimized  if  At  <  1/bandwidth.  Since  Uj(f)  is  repetitive  or. 


55 


the  interval  f, 


it  is  useful  to  define  the  time  function 


*  • 


corresponding  to  only  one  period: 

A 

J9(f  )expj2Hft  df 


u(T)  =  j  Ui 


(77) 


This  leads  to  the  relationship: 


and 


■/ 


u(kAt  )  =  /  ( f )expj2nfkAt 


4  00 


us(f)  =£^Atu(kAt  )exp-  jZTlfkAt 


(78) 


(79) 


For  any  finite  duration  sequence,  the  discrete  Fourier 
transform,  X(K),  of  the  sequence  x(n)  is: 


X  (  K )  = 


A/*  I 


»*o 


x  ( n )  exp(  -  j  2Ttkn/N  )  ,  0<k<N-1 

0  ,  e.w. 


(80) 


N  is  the  number  of  samples  in  the  sequence.  The  inverse 
transform  is: 

1 /Ny~X(K)exp( j2Hkn/N) ,  0<n<N-1  (81) 

*»©  o  ,  e.w. 

Equations  81  and  82  can  be  implemented  via  computationally 
efficient  algorithms  such  as  Fast-Fourier  Transforms  ( FFT ) . 
Also,  discete  convolutions  ( ie  z ( t ) =y ( t ) *x( t ) )  can  be 
performed  by  computing  the  Y(K)  and  X(K)  sequences,  computing 
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IV  Computer  Simulation  Development 


Approach 

The  objectives  of  this  study  are  to  generate  the 
following  for  each  signal  tested: 

(1)  Ambiguity  diagram  plot 

(2)  Signal  spectrum  plot 

(3)  Received  signal  spectrum  plot  (due  to  10 
scatterer  target) 

(4)  Ideal  matched  filter  output  plot 

(5)  Target  filter  spectral  characteristics  plot 

(6)  Incoherent  target  filter  output  plot 

(7)  Coherent  target  filter  output  plot 


Objectives  one  through  four  are  included  to  verify  the  work  of 
Salzman  (Refll)  and  provide  a  complete  description  of  the 
different  signal  processing  and  target  modeling  schemes 
simulated.  Objectives  six  and  seven  are  the  main  results  which 
this  paper  requires.  The  data  generated  for  these  two  plots 
will  be  compared  at  the  appropriate  sample  time  in  order  to 
determine  which  filter  implementation  provides  superior 
signal- to-noise  ratio. 

In  this  section  of  the  paper,  the  algorithms  developed 
will  be  explained.  Appendix  A  contains  the  system  flowchart 
for  the  software  and  appendix  B  contains  the  program  listing. 
The  software  is  written  in  FORTRAN  V. 


Calculat ing  the  Magnitude  of  the  Ambigu ity  Function 

The  ambiguity  function  (and  diagram)  is  the  correlation 
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of  the  complex  envelope  of  the  signal  with  a  doppler-shifted, 
time-delayed  version  of  itself;  that  is,  the  ambiguity 
function  is  the  two-dimensional  correlation  function  in  delay 
and  doppler.  It  is  defined: 

'y.  0?,  w )  =  J"u(  t  )u  *(  t->)  exp[j2fivt]  dt  (83) 

-  cm 

The  ambiguity  function  can  be  expressed  as  a  convolution  of 
the  following: 

y  ( r,  w )  =  u(>)  exp[  j  2ttv>]  *  u*(-?1  (84) 

In  this  simulation  we  define  the  arguments  of  the  convolution 
as  finite  sequences  and  calculate  the  ambiguity  function  via  a 
linear  convolution.  Each  argument  is  sampled  to  form  an  N- 
point  sequence  and  then  at  least  N-1  zeros  are  added.  The  N-1 
zeros  are  necessary  for  the  discrete  convolution  calculations 
to  be  correct  (Ref 4 : 85 -88 ) .  The  linear  convolution  is 
calculated  by  the  indirect  method  which,  given  two  time 
sequences  x#  and  y*  ,  consists  of: 

1  )  compute  the  sequences  X^  and  Y*  via  the  discrete  Fourier 
transformation  ( FFT )  of  xK  and  y(<  respectively. 

2)  form  Z^  =  X^Y^  for  X=0,1,2,...2N-1 

3)  compute  the  sequence  z «  via  the  discrete  Fourier 
transform  (FFT)  of  Z*  (Ref4:98) 


Sufficient  sampling  of  the  waveform  requires  knowledge 


of  the  Nyquist  rate.  The  time-bandwidth  product  (TB  or  signal 
duration  multiplied  by  signal  bandwidth)  of  the  complex 
envelope  provides  the  Nyquist  rate  (Ref6:254).  If  the  signal 
is  sampled  at  1/TB  or  less,  the  N-point  sequence  is  generated. 
The  Fast  Fourier  Transform  algorithm  used  in  this  simulation 
was  available  through  the  IMSL  mathematical  library  (CDC  Cyber 
computer  system).  Good  resolution  of  the  resultant  convolution 
may  require  sampling  rates  which  are  higher  than  the  Nyquist 
rate.  The  sequences  generated  by  the  software  are  limited  to 
1024  points.  This  constraint  was  primarily  imposed  by  the 
100,000  words  of  compilation  memory  available  for  interactive 
computations  on  the  Cyber.  This  limitation  will  severely 
constrain  both  the  number  of  pulses  and  subpulses  which  can  be 
simulated.  This  in  turn  limits  the  length  of  coding  (either  FM 
or  phase)  possible  and  will  be  elaborated  on  when  a  simulation 
waveform  is  selected. 

Sampling  repeated  and/or  coded  waveforms  requires 
careful  consideration.  Each  element  of  the  code  ( ie  subpulse 
interval)  must  be  sampled  equally  in  order  to  eliminate  any 
bias  and  skewing  of  the  plotted  function.  Consider  a  two  pulse 
signal  with  each  pulse  1  second  in  duration  and  duty  factor 
of  1/3.  If  we  choose  to  make  16  samples  of  the  waveform,  each 
sample  is  separated  by  3/8  seconds.  This  results  in  two 
samples  made  of  the  first  pulse  and  three  samples  made  of  the 
second  pulse.  This  would  of  course  result  in  an  incorrect 
answer.  The  problem  is  avoided  by  weighting  the  samples  in 
such  a  manner  as  to  insure  each  elemental  portion  of  the 
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waveform  is  sampled  uniformly.  We  have  specified  the  number  of 
samples  we  wish  to  make,  the  signal  duration,  and  pulse  width. 
By  forming  the  expression  [number  of  samples/signal  durationjx 
[pulse  width]  and  truncating  the  result,  we  obtain  the  largest 
integer  value  for  the  number  of  samples  per  pulse  which  can  be 
achieved  and  insure  uniform  sampling  (Ref8:45). 

As  an  example,  assume  we  are  simulating  a  phase  coded 
waveform  where  the  amplitude  of  the  complex  envelope  is 
modulated.  Given  a  single,  1  second  pulse  of  unity  amplitude, 
we  choose  to  make  8  samples.  Also,  assuming  a  5-element  Barker 
code,  the  phase  of  each  element  would  be  0,0,0, 1^,0.  The 
sequence  corresponding  to  u(T)exp[jw"M  would  be: 

exp[ jw( 1 / 5 ) ] ,exp[ jw{ 2/5) ] ,exp[ jw( 3/5 )] , -exp[ jw( 4/5 ) ] 
exp[ jw] ,0,0, 0,0, 0,0, 0,0, 0,0,0. 

The  sequence  representing  u*  (  ~"V)  is: 

0,0, 0,0, 0,0, 0,0, 0,0, 0,1  ,-1 ,1 ,1  ,1 


Once  the  proper  sequences  are  formed  for  u(T)exp[jwJ"  ) 
and  u  ( ,  the  linear  convolution  must  be  performed.  This  is 
accomplished  via  the  indirect  method  as  previously  stated.  The 
simulation  calculates  the  ambiguity  diagram  values  as  a 
function  of  delay  given  a  particular  doppler  value.  The  result 
is  a  cross-sectional  view  of  the  diagram  which  is 


perpendicular  to  the  doppler  axis.  Any  value  of  doppler  may  be 


selected,  however,  this  study  will  always  select  the  zero 
doppler  cut  since  this  corresponds  to  the  output  of  a  matched 
filter  which  is  perfectly  matched  to  the  waveform.  In  this 
manner,  we  are  only  viewing  a  particular  cross-sectional 
"slice"  of  the  total  ambiguity  diagram.  After  calculating  the 
values  of  the  ambiguity  function,  their  magnitude  (or 
magnitude  squared)  is  plotted.  This  plot  represents  the  best 
possible  range  response  of  the  radar  waveform. 

All  of  the  ensuing  sections  of  this  chapter  describe 
software  which  rely  extensively  on  the  algorithm  just 
described.  The  same  proceedure  for  setting  up  the  proper 
sampling  and  resultant  sequences  is  fundamental  to  matched 
filter  and  target  filter  calculations  as  well.  This  will 
become  apparent  as  the  chapter  develops. 

Signal  Spectrum  Calculations 

Two  signal  spectrum  calculations  (plots)  are  generated. 
The  first  is  simply  the  signal  spectrum  of  the  original  radar 
waveform,  and  the  second  is  the  spectral  characteristics  of 
the  received  signal.  The  first  is  generated  by  performing  a 
discrete  Fourier  transform  of  the  sequence  corresponding  to 
u( ^ )  and  plotting  the  result.  Generation  of  the  received 
signal  spectrum  is  not  as  straight  forward  and  requires  the 
following  assumptions:  1)  the  target  consists  of  ten 
independent  scatterers  which  may  or  may  not  be  fluctuating  in 
amplitude,  2)  the  target  acts  as  a  filter  and  thus  can  be 
described  via  a  range  extended  scatterering  function,  3)  the 


range  extended  scatterering  function  can  be  modeled  as  a 
decaying  exponential  in  the  frequency  domain  (Ref5).  The 
procedure  is  to  specify  (in  frequency)  the  exponential  as  a 
Gaussian  pulse  centered  at  zero.  From  this  model  of  the  range 
extended  scatterering  function,  a  complex  impulse  response  is 
developed.  The  original  signal  spectrum  is  multiplied  by  the 
complex  range  extended  scatterering  function,  Fourier 
transformed,  and  plotted.  The  development  of  the  range 
extended  function  starts  with  a  real,  bandpass  Gaussian 
spectrum: 


H(jw)=  exp[  -  tl(  w-we  )a  ]  +  exp[ -*(  w+w0  )Z  ] 


(85) 


noting  that, 


/V*  . 


H(jw)=  H(w-w0)  +  H  ( -w-wc  ) 


(86) 


where  H(jw)  is  complex  and 

H(jw)=  5f[h(t)]  =  Re  [h(  t  )exp(  jwD  t )  ] 

A#  rs* 

=  h(  t )  exp(  jwo  t )  +  h  ( t )  exp(  -  jw„  t ) 


(87) 


Equation  85  dictates  the  real  impulse  response: 


h(t)=  1  / 2jTf3'  exp[-tl/4o(  cosw*  t 


thus, 


H( jw ) =  1/2 


1  / 2>jTT«(  exp[ -t*  /4*]  exp[  j  w©  t  j  + 


1/2  77*  exp[ -t'  / 4 -<] exp[  -  jwe t.  ]  (88) 


and  by  comparing  equations  88  and  87  we  see  that 


h(t)  =  1/2JW3  exp[  -t*  /  4«l ) 

(89) 

which  implies 

that 

H(jw)=  exp[-*w2] 

(90) 

Equation 

89  requires  the  addition  of  a  time  shift. 

t0 ,  which 

accounts 

for 

causality  (Ref2:8c),  and  results  in 

the  final 

form  of  the  complex  range  extended  scatterering  function: 

H(jw)=  exp[ -o(wa  ]exp[ -jwt0  ]  (91) 

Equation  91  is  simulated  directly.  $(jw)  is  sampled  according 
to  the  same  criteria  as  the  signal,  multiplied  sample  by 
sample  with  the  original  signal  frequency  domain  samples,  and 
plotted  as  the  received  signal  spectrum.  The  process  is 
repeated  for  both  fluctuating  and  non-fluctuating  target 
cross-section. 

Matched  Filter  and  Target  Filter  Calculations 

The  procedure  for  calculating  both  the  target  and 
matched  filter  outputs  is  essentially  the  same.  The  only 
difference  is  that  a  unique  filter  impulse  response  is 
simulated  for  each. 

The  matched  filter  output  is  the  discrete  convolution 
of  the  received  signal  time  samples  with  a  conjugated,  time 
inverted  version  of  the  original  signal.  The  next  section 
develops  the  governing  equations. 


The  target  filter  output  is  the  discrete  convolution  of 
the  received  signal  with  the  Fourier  transform  of  the 
conjugated  range  extended  scatterering  function  impulse 
response.  The  target  filter  impulse  response  parameter  ol  , 
generally  is  not  equal  to  the  d  which  characterizes  the 
received  signal  structure.  Specifically,  equation  89  or 
alternatively  91,  is  used  to  generate  the  received  signal  ( ie 
model  the  target)  and  designate  the  target  filter  impulse 
response, but  different  values  of  ol  are  used  in  each.  A  larger 
value  of  ol  will  be  used  for  the  target  filter  impulse 
response,  resulting  in  a  narrowbanding  of  the  return.  This  is 
a  critical  point  in  that  narrowbanding  the  return  may  provide 
and  is  indeed  necessary  to  yield  a  practical  spread  spectrum 
system  (assuming  current  technology  limitations  as  presented 
in  chapter  II ) . 

As  a  last  note,  it  is  apparent  that  we  have  assumed  a 
received  signal  which  was  generated  via  a  known  filtering 
process  ( ie  the  range  extended  scatterering  function  of 
equation  91),  and  further  that  we  have  set  up  a  matched  filter 
(the  target  filter)  which  compensates  for  the  known  target 
scatterering  characteristics.  This  might  appear  to  be  a  self- 
fulfilling  prophesy  or  a  circular  argument.  Two  important 
points  must  be  considered:  first,  the  validity  of  the  results 
rests  partly  on  the  validity  of  the  assumption  and 
approximation  of  the  form  of  the  ten  scatterer  target  density 
function.  For  this  assumption,  the  author  has  relied  on  the 
article  "Effect  of  a  Few  Dominant  Specular  Reflectors  Target 


Model  Upon  Target  Detection"  (Ref 5:670-671 ) .  Second,  the  value 
of  for  the  target  filter  is  larger  than  the  one  used  for 
generation  of  the  received  signal,  thus  constituting  the 
ability  or  inability  to  successfully  narrowband  the  return 
prior  to  range  and  doppler  processing.  With  the  above  in  mind, 
the  next  section  will  develop  the  equations  which  support  this 
section  and  define  the  relationship  between  coherent  and 
incoherent  processing. 


Coherent  and  Incoherent  Target  Filtering 

This  section  will  develop  the  equations  for  coherent 
and  incoherent  filtering  of  complex  signal  envelopes.  It 
should  be  kept  in  mind  that  the  target  filter  is  actually  a 
matched  filter  which  is  matched  to  the  received  signal. 

On  page  338  of  DiFranco  and  Rubin  (Ref3)  we  find: 
"Matched  filter  processing  for  a  pulse  train  is  sometimes 
called  coherent  processing  or  predetection  integration".  In 
other  words,  the  matched  filter  is  inherently  a  coherent 
device.  The  phase  matching  is  based  on  the  signal  spectrum 
alone  and  results  in  a  maximum  peak  response  occuring  at  t=T, 


where  T  is  the  causality  delay.  The  amplitude  matching  is 
dependent  on  both  the  signal  and  noise  spectrum,  resulting  in 
the  maximization  of  the  signal-to-noise  ratio  (Ref 2:On 
Understanding  the  Matched  Filter  in  the  Frequency  Domain  by 
Theodore  G.  Birdsall).  The  complex  response  function  is 


*♦  A© 


(>(**)=  1/2  J  *<t)  vk*(t-»  dt 


(  92  ) 
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where  4>(t)  =  s(t)+js(t) 
and  s(t)  is  the  real  signal. 


s(t)  is  the  Hilbert  transform  of  the  real  signal. 


Further, 


*  '  -m* 

He[f(>)]=  1/2  /  s(t)s(t-T)  dt  +  1/2  /  s(t)&(t-M  dt 

.  J  -*>• 


(93) 


In  real  notation,  the  real  autocorrelation  function  (matched 
filter  response)  of  s(t)  is 


r(» 


s(t)s(t-?*)  dt 


■/ 


S(f)|  exp[j2TfM  df  (94) 


s(t)  has  the  same  spectrum  as  s(t),  thus  the  integrals  of 
equation  93  are  equivalent  (Ref9:32). 

We  note  that 

r( ?■)  =  Re  [  p(7)] 

Also, 


4>(t)  =  u(  t  )exp[  j21tf#?’] 
where  fo  is  the  center  frequency  and 


f  (» 


1/2 


exp[  j  21»f0 


J u(  t)u*  (t-»  dt 


(95) 


Equation  95  is  simulated  directly  (convolution  process 
discussed  previously)  with  f0  =0  by  ignoring  the  insignificant 
conversion  to  video.  The  actual  plotted  response  is  calculated 
by  taking  the  real  part  of  equation  95. 


The  incoherent  target  filter  is  the  combination  of  the 
matched  filter  followed  by  an  envelope  detector (Ref 1 2: 201 ) .  In 
other  words,  the  phase  information  is  destroyed  and  the 
samples  are  incoherently  collapsed  (added).  Starting  with 
equation  95,  the  envelope  of  the  real  target  filter  response 
(incoherent)  is: 

Env |  r(?-)^=  |f{?-) 

-  1/J/u  (t)u’(t-r)  dt  (96) 

where  again  exp[jw0t]  has  been  suppressed.  Equation  96  is 
simulated  via  the  techniques  of  the  previous  sections  of  this 
chapter. 

In  summary,  both  target  filters  are  simulated  by 
discrete  convolutions  according  to  equation  95  with  f©  =0  and 
u*(  t-t)  replaced  with  the  range  extended  scatterering  function 
impulse  response.  The  coherent  target  filter  takes  the  real 
part  of  this  result,  where  as  the  incoherent  filter  requires 
taking  the  magnitude  of  the  integral.  These  two  output  are 
compared  at  the  sample  point  corresponding  to  the  V=0,  v=0 
point  on  the  ambiguity  surface.  The  larger  value  constitutes 
the  greater  signal-to-noise  ratio  at  the  output  of  the  since 
r(0)=  p(0)=energy  (Ref9:33).  The  next  chapter  presents  the 
results  of  the  processes  presented  above. 


V.  Software  Application  and  Results 


Simulaton  Test  Waveform 

Due  to  computer  memory  limitations  and  thus  array  size 
constraints,  the  maximum  number  of  data  points  is  fixed  to 
1024.  With  this  in  mind,  a  four  pulse,  rectangular  pulse  train 
of  1/4  duty  factor  was  selected  as  the  base  signal.  The  signal 
is  then  coded  in  phase  or  frequency  as  will  be  discused.  The 
above  selection  results  in  a  maximum  of  512  samples  (after 
weighting)  of  the  signal  and  up  to  32  allowable  subpulses  or 
elements. 

The  key  to  the  simulation  process  is  the  ambiguity 
calculations  and  resultant  plot.  All  other  matched  filter, 
signal  spectrum,  and  target  filter  plots  follow  the  same  basic 
software  manipulations.  A  single  run  through  the  algorithm 
produces  one  of  the  following: 

1)  Ambiguity  plot  and  signal  spectrum  plot. 

2)  Received  signal  plot  and  matched  filter  plot. 

3)  Target  filter( coherent  or  incoherent)  spectrum 
plot  and  target  filter  ouput  plot. 

In  order  to  test  the  software,  a  13-element  Barker  code  was 

injected  into  the  simulation.  The  verification  signal 

contained  13  pulses,  phase-coded,  with  1/2  duty  ratio.  The 

resultant  spectrum  and  ambiguity  cut  (doppler=0)  are  shown  in 

figures  13  and  14  respectively.  Both  plots  are  verified  in  the 


literature , 


in  particular  by 


Sa Izman ( Ref  1 1  ) 


Rihaczek (Ref9:215-21  6)  . 


Phase  Coded  Signal 


The  basic  test  waveform  of  the  previous  section  was 
intrapulse  phase  coded  (0,7T)  and  simulated.  Each  of  the  four 
pulses  contain  25  elements  generated  from  a  pseudo-random  code 
of  length  63.  A  50  element  section  of  the  code  is  repeated 
throughout  the  signal.  Figure  15  is  the  signal  spectrum  and 
figure  16  is  the  zero-doppler  cut  of  the  phase  coded  signal. 
The  central  peak  of  figure  16  designates  the  sample  point  of 
the  matched  filter  and  is  located  at  15.96  on  the  time  units 
scale.  Note  the  level  of  the  two  ambiguious  peaks  on  either 
side  of  the  main  response  peak.  This  provides  a  clear  area  of 
approximately  14  time  units  (from  9.0  to  23.0).  The  clear  area 
is  that  delay  span  in  which  signals  can  be  searched  and 
separated  in  range  ( ie  area  of  sufficiently  low  ambiguious 
peaks).  The  central  peak  is  at  40dB,  with  the  next  highest 
response  level  within  the  clear  area  at  approximately  21 dB. 
This  is  clearly  a  thumbtack-type  ambiguity  function  as  would 
be  expected  from  a  coded  or  noise-type  radar  waveform  (ReflO). 
Figures  17  and  18  show  the  received  signal  spectrum  and 
matched  filter  output  for  non-fluctuating  target  cross- 
section.  The  sample  point  has  a  magnitude  of  15.59dB.  Figures 
19  and  20  show  the  target  filter  spectral  characteristics  and 
coherent  target  filter  outputs  respectively.  The  coherent 
target  filter  has  been  modeled  with  the  same  value  of  d  as  the 
received  signal  spectrum.  The  magnitude  at  the  sample  point 
(15.96  time  units)  is  7.58dB.  The  corresponding  value  for  the 
incoherent  target  filter  was  9.34dB.  Thus,  the  incoherent 


filter  shows  a  1.76dB  improvement  in  output  signal-to-noise 
ratio  for  the  non-f luctating  case. 


The  following  series  of  plots  are  for  the  coherent 
filter(s).  The  incoherent  filter  output  is  identical  except 
for  a  slight  widening  of  the  main  lobe.  Figures  21  and  22 
depict  the  received  signal  spectrum  and  matched  filter  output 
for  scan-to-scan  fluctuating  RCS  target  (Swerlingl).  The 
effect  appears  to  be  slight  and  further  RCS  fluctuations  will 
be  pules-to-pulse  or  Swerling  II  modeled.  Figures  23-26  show 
the  same  series  of  events  but  e(  has  been  reduced  from  0.1  to 
0.01.  The  coherent  target  filter  does  appear  to  widen  the  main 
lobe  response  slightly.  The  next  sequence  of  plots,  figures 
27-30,  are  for  a  target  of  different  Swerling  II  statistics. 
The  results  are  similar  to  the  previous  plots.  In  figures  31- 
34,  the  value  of  ot  used  to  generate  the  received  signal 
spectrum  is  0.01,  whereas  the  value  for  the  target  filter 
spectrum  equals  0.01.  The  result  of  narrowbanding  the  return 
by  60%  is  the  smooth  response  shown  in  figure  34.  The 
incoherent  target  filter  performs  similarly  with  a  sampled 
value  7.62dB  higher.  Figures  35-38  are  the  same  sequence  of 
events  once  again  with  o<  (received  signal  spectrum)  =  0.00' 
and  cA  (coherent  target  filter)  =  0.2.  The  results  for  these 
plots  (including  data  for  the  incoherent  filter)  and  an 
additional  20  runs  are  tabulated  in  table  II. 

It  is  clear  from  table  II  that  the  coherent  target 
filter  produces  too  narrow  a  main  peak  and  thus  an 
unacceptable  sampled  ouput  signal  power  for  values  of  'X  0.1. 
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Further,  in  each  of  the  simulation  runs  the  incoherernt 
realization  provided  the  best  output  signal- to-noise  ratio. 


Interpulse  Frequency  Hopping  Waveform 

A  four  pulse,  20%  duty  cycle  signal  was  simulated  over 
10  runs  with  Swerling  II  RCS  fluctuations.  The  results  are 
tabulated  in  table  III.  Again,  the  coherent  target  filter 
proved  slightly  superior.  The  ratio  of  incoherent  to  coherent 
signal  output  power  remained  within  O.OIdB  for  the  10  run  data 
set . 

Interpulse  Frequency  Hopping  with  Subpulse  Chirp 

The  same  basic  waveform  was  again  simulated  with 
intrapulse  chirp  (25  seperate,  linear  frequency  slopes)  added. 
The  incoherent  target  filter  once  again  shows  a  slight 
advantage  in  output  signal  power.  Table  IV  shows  the  results 
for  5  separate  runs  with  Raleigh  (Swerling  II)  RCS 
fluctuations.  The  non-fluctuating  scenario  is  also  included. 
The  constant  average  advantage  of  the  incoherent  filter  varies 
only  a  few  tenths  of  a  decibel. 

Summary 

It  appears  that  for  all  three  waveforms,  and  under 
Swerling  II-type  RCS  fluctuations,  the  incoherent  target 
filter  is  superior  in  terms  of  greater  output  signal  power. 
This  will  yield  a  higher  ouput  signal - to-noise  ratio.  The 
advantage  of  the  incoherent  filter  is  generally  on  the  order 


Target  Filter  Sampled  Value(dB) 


Coherent  J 

Incoherent 

★ 

RCS  fluctuation 

0(  (filter) 

\ 

7.58 

9.34 

1  .76 

non -fluctuating 

- 

-60.00 

-0.32 

59.70 

Swerling  II 

0.01 

-60.00 

0.44 

60.40 

Swerling  II 

0.01 

-7.87 

-0.25 

7.40 

Swerling  II 

0.10 

-2.92 

1  .00 

3.90 

Swerling  II 

0.20 

-2.84 

-1  .32 

1  .  50 

Swerling  II 

0.20 

-1  .83 

-0.47 

1  .30 

Swerling  II 

0.20 

-1 .85 

0.09 

1  .90 

Swerling  II 

0.20 

-3.75 

-1  .95 

1  .  80 

Swerling  II 

0.20 

-3.52 

0.24 

3.80 

Swerling  II 

0.20 

0.27 

1.15 

0.90 

Swerling  II 

0.30 

-5.36 

-3.63 

1  .70 

Swerling  II 

0.30 

0.67 

0.97 

0.30 

Swerling  II 

0.40 

-3.92 

-2.39 

1  .50 

Swerling  II 

0.40 

-2.18 

-1  .66 

0.50 

Swerling  II 

0.40 

-2.42 

-1  .45 

0.90 

Swerling  II 

0.40 

-0.79 

-0.24 

0.50 

Swerling  II 

0.  10 

-1  .62 

-0.98 

0.60 

Swerling  II 

0.40 

-0.57 

-0.18 

0.40 

Swerling  II 

1  0.40 

-0.58 

-0.12 

0.50 

Swerling  II 

0.40 

-2.61 

-2.00 

0.60 

Swerling  II 

|  0.40 

-0.22 

0.43 

0.70 

Swerling  II 

|  0.40 

-1.46 

-1.00 

0.50 

Swerling  II 

,  0.40 

-3.25 

-2.37 

0.90 

,  Swerling  II 

0.40 

-5.  35 

-4.49 

0.90 

j  Swerling  II 

0.40 

*  -  the  ratio  of  incoherent/coherent  target 
filter  output  signal  power  in  dB. 


Table  II.  Coherent  vs.  Incoherent  Target  Filter  Output- 
Signal  Power  for  a  Phase  Coded  Waveform 


Coherent 


-5.98 

-4.41 

-9.78 

-11.49 

-5.60 

-5.03 

-10.20 


Incoherent 


-6.64 

-5.40 

-6.60 

-5.24 

-3.67 

-9.06 

-10.76 

-4.85 

-4.28 

-9.45 


RCS  Fluctuation 


4  non- fluctuating 

5  Swerling  II 

5  Swerling  II 

0.74  Swerling  II 

0.74  Swerling  II 

0.72  Swerling  II 

0.73  Swerling  II 

0.75  Swerling  II 

0.75  Swerling  II 

0.75  Swerling  II 


(filter ) 


0.40 

0.40 

0.40 

0.40 

0.40 

0.40 

0.40 

0.40 

0.40 


*  -  the  ratio  of  incoherent/coherent  target 
filter  output  signal  power  in  dB. 

Table  III.  Coherent  vs.  Incoherent  Target  Output  Signal 
Power  for  Interpulse  FM  Waveform 


Target  Filter  Sampled  Value  (dB) 


Coherent 

Incoherent 

19.24 

19.42 

14.62 

14.18 

16.52 

16.70 

14.45 

14.63 

15.97 

16.15 

15.83 

16.02 

RCS  Fluctuation 

non-f luctuating 

Swerling 

II 

Swerling 

II 

Swerling 

II 

Swerling 

II 

Swerling 

II 

( filter ) 


the  ratio  of  incoherent/coherent  target  filter 
output  signal  power  in  dB 


Table  IV.  Coherent  vs.  Incoherent  Target  Filter  Output 
Signal  Power  for  Interpulse  w/subpulse  Chirp 
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VI  Conclusions  and  Recommendations 

Conclusions 

The  objective  of  this  study  was  to  determine  whether  a 
coherent  or  incorerent  target  filter  provided  greater  output 
signal-to-noise  ratio.  The  signals  under  test  were  all  coded, 
spread  spectrum  waveforms  designed  for  use  in  a  coherent 
airborne  radar  system.  It  was  noted  that  the  computer 
simulation  was  dependent  on  the  characterization  of  the  range 
extended  scatterering  function.  Knowing  the  scatterering 
characteristics  and  the  accuracy  of  their  approximatons  ( ie 
decaying  exponential,  Gaussian  pulse)  is  paramount  to  the  type 
of  analysis  performed. 

The  25  element,  pseudo-random  pase  coded  signal 
analysis  showed  two  main  features.  First,  values  of  o(  less 
than  0.1  yielded  50-60  dB  higher  ouput  signal  power  for 
incoherent  processing.  Values  of  U  at  or  above  0.1  allowed  a 
smoother,  wider  peak  around  the  15.96  time  unit  sample  point. 
Higher  values  of  o t  correspond  to  a  narrower  filter  bandwidth. 
The  consistent  performance  edge  (1-3  dB)  of  the  incoherent 
target  filter  found  by  this  study  agrees  with  the  findings  of 
Salzman  (Ref 1 1 :69-70) .  Second,  the  incoherent  target  filter 
provides  improved  signal-to-noise  ratio  due  to  its  broader 
main  peak.  This  produces  a  higher  sampled  value  away  from  the 
output  centroid  (located  at  15.96  time  units)  as  compared  with 
the  coherent  filter. 

The  interpulse  frequency  hopping  and  frequency  hopping 


with  subpulse  chirp  show  similiar  results.  In  both  cases,  the 
incoherent  filter  is  superior.  The  performance  difference  for 
these  two  waveforms  is  much  more  consistent  as  compared  to  the 
phase  coded  signal.  Interpulse  coding  yields  a  steady  0.72- 
0.75  dB  incoherent  filter  advantage,  whereas  interpulse  FH 
with  subpulse  chirp  holds  around  0.4 -0.7  dB  improvement 
margin. 

It  appears  that  with  assumptions  made  concerning  RCS 
fluctuations,  target  scatterering  characteristics,  and  radar 
waveform,  that  the  incoherent  target  filter  is  superior  to  the 
coherent  target  filter  on  the  basis  of  output  SNR. 

Recommendat ions 

In  addition  to  the  signals  tested,  the  simulation 
software  is  capable  of  testing  linear  frequency  modulated, 
variable  frequency  modulated,  intrapulse  frequnecy  hopping, 
and  intrapulse  FH  with  subpulse  chirp  waveforms.  The  software 
is  also  capable  of  simulating  Gaussian  as  well  as  rectangular 
pulse  trains.  Further,  different  types  of  RCS  fluctuations  and 
target  scatterering  characteristics  can  be  easily 
accommodated.  Due  to  time  constraints  involved  in  producing 
this  paper,  these  possibilities  were  not  explored.  Whether  or 
not  they  provide  further  insight  into  the  applicability  of 
the  target  filter  is  left  for  further  reasearch.  However, 
based  on  the  data  presented,  an  incoherent  target  filter  would 
be  advantageous  when  SNR  is  a  factor  and  ease  of 
implementation  important. 


Future  studies  might  consider  testing  those  waveforms 
not  tested  in  this  study  with  and  without  different  range 
extended  scatterering  functions.  Alternatively,  the  concept  of 
non-cooperative  target  recognition  (NCTR)  might  be  addressed 
by  extending  the  software.  NCTR  is  an  attempt  to  exploit  the 
range  super-resolution  phenomenon  (rather  than  compensate  for 
it)  in  order  to  provide  a  type  of  target  signature.  In  this 
manner,  the  target  could  be  detected,  identified,  and  classed. 
The  software  used  in  this  study  might  be  extended  to 
investigate  target  super-resolved  range  signatures. 
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Ppectnar  Analysis,  Matched  Filter,  and  Tarpet 
Filter  Processing 


PJ.r  Progran  'MAIN* 


cal  1  ’Input’  to  defir.e 
pulse  parameters 


call  ’Code’  to  set 
up  pseudo-ebde  values 


/  13  X 

'pulse  phase  \ 
shift  coded 


call  ’Phase' to  Ret 
frequency  and  .arpHtui; 
for  vaefonn  a.molcu 


ceTVfM'  to  pat  frequency 
:  u  la  tier,  of  waveform 


>11  ’Recpl’  to  yet 
amplitude  of  ".veform 
Miupl  i 


call  ’Setup’  to  pit 

:;e  uence  ;  in  proper  t\,rr-o : 
for  11  -ear  cor.v. letters 


Copy  available  to  DTIC  does-  tv 
oetmii  felly  legible  r?p*edt!  ‘it:. 


a  •clrul  t  Veii.g 


gene-r  tel 


PROGRAM  MAIN 


100= 

110= 

120= 

130=C  PROGRAM  MAIN 

1 40=C  WRITTEN  BY  KENNETH  W.  ALBERT/ < REF l Rt ED) 

150=C  17  OCT  03 

160= 

170= 

1 80=  COMPLEX  XVALS  <1024 ) .  YVALS < 1024 ), Z  <1024  > ,  VAI  UE ,  lih  LAY 

190=  INTEGER  SAMPLE, PERIOD! 1024) ,D,C<64 ) 

200=  REAL  MAG< 1024) ,FREQ! 1024) , FILTER! 1024) .MMAG! 1024) ,PRT( 50) 

210=  OPEN! 10,FILE= 'DATA' ) 

220-  OPEN<  12,FIL.E='RESni  2' ) 

210=  OPEN!  1 5, FILE=' RESULT  > 

240=  OPEN <11. FILE* 'RFSULi'  ) 

250-  REWIND  10 

260-  RFWIND  11 

270=  REWIND  12 

280-  REWIND  15 

290= 

300= 

310=  CALL  INPUT < PERIOD , PRT , PW , NPRT , TOT Al. , SAMPLE , NUM „  M ,  N ) 

320=  CALL  CODE ( C ) 

330  RfcAIKlO,*)J 

340=  IF  (J.EQ.l)  GO  fO  5 


350  (AIL  PHASE !FRFQ,PRT  .PERIOD,?  ,MfN,N-'RT  ,GFACT  , 

360  HPW, SAMPLE, TUT Al.  ,C> 

370  bU  TO  20 

3R0-5  URTTE<*r«>/Pm.SE(S/  IS  FREQUENCY  CODED 

390  -  CALL  FiKFREO, PERIOD. NPKT.PW,N,C,M,TOIA!  ,PRT> 

400  CALL  RECPL (PERIOD, 7, br ACT, NPRT,N> 

410-  20  SCAI  E=C»FACT/!N*NPRI ) 

420=  WRITE! 12, *)M 

430=  WRITE!15,*)M 

440=  WRITE! 11, *)M 

450  -  RFAD<10,*)W 

460-25  VALUE=CMPLX!O.C,W) 

470=  WRITF!*,300)W 

480=  READ! 10,4 )D 

490=  WRITE!*,*) 'THE  VALUE  OF  D  IS',D 

500=  IWK=6*SAMPLE+150 

510=  WK=6*SAMPIF  P150 

520=  CALL  SETUP ! 0 , Z , FREQ , XVALS , YVAL8 , OAI  UP , : U I AL • ->  AfiPL I 

5  30=  BIWK,  WK,PW,fl) 


540= 

CALL  FFTCC ! YVALS , SAMPLE , IWK , WK  ) 

550- 

IF(D.EQ.0)THEN 

560- 

WRITE! *,301) 

570- 

ELSE 

580- 

TF!D,EQ.l)  THEN 

:  90  = 

WRITE!*, 307) 

600= 

ELSE 

•5  0  = 

WRITE!*. 315) 

a^v  9 


P ftOrtt 


620= 

FND  IF 

630= 

END  IF 

640= 

DO  30  1=1, 2*M 

650= 

TIME=< TOTAL/M)*! 1-1) 

660= 

X=CABS!YVALS! I>> 

670= 

IFGX.LT. 0.001)  THEN 

680= 

X=0»00l 

690= 

END  IF 

700= 

MMAG (I )  =  1 0* ALOG 1 0  ( X *X ) 

710= 

IF(D.EG.0)THEN 

730= 

WRITE! 15,305>TIME  »MMAG<  I ) 

740= 

END  IF 

750= 

IF!D.EQ.1>  THEN 

770= 

WRITE! 15,310) TIME  ,MMAG(  .1 ) 

730= 

END  IF 

790=30 

CONTINUE 

810= 

IF  <D.  EG).  3)  THEN 

815= 

READ! 10,*) ALPHA? 

820= 

DO  42  1=1, SAMPLE 

830= 

TIME=(PW/N)Kt(I-l> 

840= 

DELAY=CMPLX!0 .0,62. ) 

050= 

FILTER! I )=EXP!-Al  PHA?*T  THF*TIMF) 

860= 

Z( I )=FILTER( I >*YVALS( I ) 

970= 

XVALS ( I ) =CON J6 ( Z ( I ) ) *CF XP !  DEL  AY *  T I HE ) 

880= 

X=CABS ( XVALS ( I ) ) 

890= 

IF (X.LT. 0.001 ) THEN 

900= 

X=0.001 

910= 

END  IF 

917= 

MMAG<I)=10*AL0G10<X*X) 

916= 

WRITE! 11, 320) TIME, AHAG< I * 

920=42 

CONTINUE 

930= 

El  SE 

940= 

CALL  FFTCC! XVALS, SAMP! E , IUK ,WK ) 

950- 

END  IF 

960= 

DO  50  1=1, SAMPLE 

970= 

XVALS ! I ) =XVAl  S ! I >  *Y VALS  < I ) 

980=50 

CONTINUE 

990= 

DO  51  1=1, SAMPLE 

1000= 

XVALS ! I ) =CONJG ! XVALS ! I ) ) 

1010=51 

CONTINUE 

1020= 

CALL  FFTCC ( XVALS , SAMPLE , I WK , UK ) 

1030= 

DO  52  1=1, SAMPLE 

1040= 

XVALS ( I ) =CONJG ! XVALS ! I ) ) /SAMPLE 

1 050=5? 

CONTINUE 

1060= 

IF!D.EQ.0)THEN 

1070= 

DO  55  1=1, SAMPLE 

1030= 

ZDUH=CABS  ( XVALS !  I ) : 

1090= 

IF(ZDUM.IT .0.001 ) THEN 

1100= 

ZD!JM=0 .001 

1  110= 

END  IF 

1 1  20= 

MAO ! 1 ) = l OKALQGl 0 <  ZDUM4ZDUM ) 

1130=55 

CONTINUE 

1140= 

ELSE 

1 1 50= 

DO  65  1=1, SAMPLE 

Hi 


1160=  Y=REAL!XVALS! I) ) 

I  170-  IK ( Y  .LT .0.001 ) THEN 

1 180-  r^O.001 

,v  1190=  END  IF 

1200=  MAGL1)=10*AL0G10!Y  * 

1  210=65  CONTINUE 

1220=  KNIT  IF 

1230=  CALI  OUT! MAG, NUM,FW, TOTAL, M,D) 

1240=  READ! 10, *>U 

1250=  IKvUI.EQ.99'?)  60  TO  500 

1 260=  00  10  25 

1270=300  FORMAT! IX. 'U-' ,F12.2) 

1280=301  FORMAT  (5X.  'TIME  DEL  A  t '  ,  5X  ,  '  b  1  6NAI  SPECTRUM  > 

I  290=305  t  DRrtAT!5X,P10.4,5X,l-  lvi.4) 

1300=307  FORMAT v FiX » ' TINE  DELAY  ,5X.  'RECIbVED  SIGNAL  SPECIfti.";  , 

1310=310  FORMAT ( 5X ,  K )  0 . 4 , 5X ,  F  J  0 . 4  > 

1320=315  FORMAT <5X, 'VALUE  OF  t ', 5X , " TARGET  FILTER  SPECTRUM'  > 
1330=320  FORMAT (5X,F10.4,5X,F10.4) 

1340=500  END 

1350=********************************************************* 
1360=**#***********#***************************#************* 

.1370=*  SUBROUTINE  INPUT 

1380=*  WRITTEN  BY  KENNETH  W.  ALBERT/ ( REF  1  REED ) 

1390=*  17  OCT  83 


1 400=***************************************yi**************** 

I  410=**************************** *********************  ******* 

1420=  SUBROUTINE  INPUT! PERIOD, PRT, PW, NPRT, TOTA;  , SAMPLE, NUN , M  ,N  * 

I  430= 

I  440=  1 N lKGE R  SAMPL  E , PK ft 1  OH ■ 1024 > 

1450=  REAl  PRT(50>  fCPRT <  50 .< 

[Jp  I  460= 

1470=  Rf nDv )0,*)PW 

1480-  READ! 10,*) I 

1490=  fK(I.EQ.O)  GO  TO  25 

1500=  WRITE! *,*) 'WAVEFORM  HAS  REPEATED  PULSES ' 

1510=  READ!  10,  *).J 

1520=  IF! J.EQ.O)  GO  TO  10 

1530=  RFAD! 10,*)SAME 

1540=  READ! 10,*)NPRT 

1550=  DO  5  1=1 ,NPRT 

1560=  PRT ( I )=SAME 


1 5/0-5  CONI  INUE 

1580-  WR I TE ! * , 300 ) SAME , NPRT 

I  590=  TOTAL =SftNf*NPRT 

1600=  00  TO  30 

In  10= 10  Wft I TE <*,*)' STAGGERED  PHIS  TO  6H  USED' 

1  ft 20=  READ! 10,*) NPRT 

1  630:-  READ!  10,*  )PRT<  1 1 

1640=  DO  15  1=2, NPRT 

1650=  READ!  !!)»*) PRT <  I ) 

1660=15  CONTINUE 

1670=  I0(AL=0.0 

1 n 80=  HO  20  l-l , NPRT 

1690=  10TAL=T0IA1  tPRl < I ) 


w 


/  •  o  •  ”  \ 


iWLHU  IV  I.Vl-*'.'’  I»V*« 


v  V'.'V  '.-  •;■  vvv 


© 


1700=  WRITE!*, 305)1, PRT!I> 

1710=30  CONTINUE 

1720=  GO  TO  30 

1730=35  WRITE! *»*) 'WAVEFORM  HAS  ONLY  0NF  PULSE' 

1740=  TOTAL =PW 

1750=  NPRT=1 

1760=  PRT ( 1 )=PW 

1770=  RFAD!10,*)NUM 

1780=  SAMPLE>2*NUh 

1790=  PERI  Oil  <  1  )=NIIM 

1800=  M=NU/1 

1810=  N=NUM 

1820=  WRITE!*, *) 'N=' ,N, 'M=' ,M, '  PERIOD! 1 )=' ,PERIGD! 1 ) 

1830=  GO  TO  36 

1840=30  READ! 10,*)NUM 

I  850=  SAMPLE=NUM*2 

1860=  IF(PW,EQ. 1 .E-ll >G0  10  32 

1870=  CPW=PW*1.F10 

1880=  CTOTAL=TOTAL*1.E10 

1890=  N=NUM/CTOTAL*CPW 

1900=  M=N*CTOTAL/CPW 

1910=  WRITE!*, *)'N=',N,'  H=',M 

1920=  DO  31  1=1 ,NPRT 

1930=  CPRT!I)=PRT(I)*1.£J0 

1940=  PERIOD! I >=N*CPRT ! I )/CPW 

1950=  WRITF!*,#)'I  =  ' ,1, '  PERIOD! I >  =  ', PERIOD! I ) 

1960=31  CONTINUE 

1970=  GO  TO  36 

1980=32  IF!PW. EG. 0.0001)  GO  TO  34 

1990=  CPW=PU* 10000.0 

2000=  CTGTAL=TOTAL* 10000,0 

2010=  N--NUM/CT0TAL*CPW 

7020-  H=N*CTOTAL/CPW 

7030=  WRITE!*, *)'N=',N,'  M=',M 

2040=  DO  33  1=1 ,NPRT 

7050=  CPRT ! I ) =PRT ! I ) *10000 . 0 

2060=  PERIOD ( I ) =N*CPRT ! I ) / CPU 

3070=  WRITE!*, *)'!='. I,'  PERIOD! 1 )=', PERIOD! I ) 

3080^;!!  CONTINUE 

2090=  GO  TO  36 

2100=34  N-NOM/TQTAL*PW 

7 110=  M=N*T0TA1/PW 

3 l 20“  WRITE!*,*) 'M=' ,M 

3130-  WR  l  T  E'  !#,*)'  N= N 

7140=  DO  35  l  l.NPRF 

2150=  PERIOD! I >=N*PRT< I )/PW 

2 1 60=  WRtlE!*,*)' 1  = ,1,  PERIOD! I >=', PERIOD! I > 

2170=35  CONTINUE 

2180=36  RETURN 

7190=300  HQRMAT! IX, 'PRTS  ARE  THE  SAME. . .PRT=' ,f J2. JO, 

3200=  @5X , ' ♦  OF  REPETITIONS*', I3> 

2210=305  FORMAT! IX, 'FOR  1= ' , I3,5X , 'PR1= ' ,F17 . 10) 

2220=  END 

7230=*****************************«*****»*«*****.t4*>.«**^f  «« 
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?240=****************************************************** 
2250=*$$$  SUBROUTINE  CODE 

22a0=****  WRITTEN  BY  KENNETH  W.  ALBERT 

2270=$$*$  17  OCT  83 

7280=****************************************************** 
2290=$$$$$$$$*$$$$$$$$$$$$$$$$$$$**$$$$$$$$*$$*$*$*$*$$$$$$$ 
2100=  SUBROUTINE  CODE<C) 

2310= 

'320=  INTEGER  X<6),C<64> 

2330 = 

2340=  X<1)=1 

7350=  DO  5  K=2»6 

7360=  X ( K ) =0 

2370=5  CONTINUE 

2380=  DO  10  1=1,64 

2390=  0< I )=X ( 1 )+X<6 ) 

2400=  IK (C< I ) .EQ.2)C( I )=0 

2410=  X ( 6 ) =X ( 5 ) 

7420=  X<5)=X<4> 

2430=  X<4)=X<3) 

2440=  X<3)=X<2) 

2450=  X ( 2)=X( 1 ) 

2460=  X(1)=C<I) 

7470=10  CONTINUE 

7480=  WRITE (*,300) 

7490=  WRITE (#,305) <C< J> » J=l ,64) 

7500=  RL TURN 

2510=300  F0RHAT(2X, ' THE  BINARY  PSEUDORANDOM  CODE  -  I.ENG1H 
7570=  \i  63  =  ') 

2530=305  FORMAT  (2X,  64  ID 
2540=  END 

2550 =*********************************************»******** 
2560=$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$*$$$$$$$$$$$$$ 

2 570=****  SUBROUTINE  PHASE 

7560=****  WRITTEN  BY  KENNETH  W.  Al  BER1/(KFK  5REED) 

2590=****  17  OCT  83 

2600=****************************************************** 

26 10=****************************************************** 
2620=  SUBROUTINE  PHASETFREQ, PRT, PERIOD, Z,M,N,NPRT,GFACT , 

2630=  (?PW , SAMPLE , TOTAL , C  > 

7640= 

>650= 

7660=  COMF1  FX  Z( 1024), AMT? 

2670=  INTEGER  PERIOD <1024), COUNT , SUM , TALLY , SAMPLE , S TAR T , F 

>680=  (*,C<64) , RATIO! 50) 

2690=  REAL  PRT < 50) , FREQ <1024) , MEAN, R <1024 ) ,CPRT < 50) 

7700=  DOUBIF  PRECISION  DSFED 

2710  - 

7720=  WRIT F<*,*) 'PULSES  ARE  PHASE-SHIFT  CODED' 

2730=  READ! 10,*)NELEM 

7740=  ITIME=N/NELEM 

2/50= ******************************************** *****>** 
7760=**  ADDITION  OF  FLUCTUATION  COEF.  MATRIX  JR 
2770=*******************************************  ****♦***»■. 


2790=  READ(10,*)F 

2790=  IF(F.EQ.l)  GO  TO  4 

'*800=  WRITE!*, *> "RAYLEIGH  AMPLITUDE  FLUCTUATIONS' 

2810=  DSEED=31313131 .DO 

2820=  DO  1  J=l, SAMPLE 

2830=  A=G6UBFS  <  DSEED  > 

2840=  B=-ALOG!A> 

2850=  R(J)=SQRT(»)/2.0 

2860=1  CONTINUE 

2870=  GO  TO  6 

2880=4  WRITE!*,*) 'NON-FLUCTUATING  CASE' 

2890=  DU  5  1=1, SAMPLE 

2900=  Fill  >  =  1.0 

2910=5  CONTINUE 

2920=*  CONTINUE 

2930:  SUMMIT I ME*NELEM 

2940  =  WRITE!*,*)' SUM= ' , SUM 

•’950--  IF(NPRT  .NE.l  )G0  TO  7 

79*0=  PERIOD!  1 »=SUM 

79  70=  WR I TF  ( * .  * »  PER  1  DTK  i  » =  ‘  ,  PER  I TJD  < 1  > 

2980=  GO  TO  30 

2990-7  IHPU.EU.l  .E-ll>bO  Hi  10 

3000=  DO  H  1=1, NPRT 

30 1 0-  l.il-;  I  •  I  » -PR  T <  1 >* l .  E 10 

3020  -  CPW=-PW*.1 .EJO 

3O30  PI  H I  illi*  1  » -SUmCPRI  i  l ) /LEW 

2040=  UP  HE!*,*)  I=',l,'  PERIOD!  1  > '  ,PF.RIUD<  I ) 

3050=8  M3NIINUE 

>060 *  Ml*  TU  30 

30/0=10  TF-PU.IVT. 0.0001)  UU  TO  20 

3080=  DO  15  I- 1, NPRT 

3090=  CPU T  U  >  =PRT ( I > *10000 . 0 

>10(1=  CPU  PW*. 1.0000,0 

3 1 10-  PERIOD! 1 )=SUM*CPRT< 1 > /CP u 

3120:  WRITE!#,  *>'!  =  ',  I,'  PERJ0D<1)=  '.PERIOD!  I ) 

3130=15  CONTINUE 

.3 1 40=  GO  TO  30 

3)50=20  DO  25  1=1, NPRT 

3  PERIOD!  I  )=SUM*PRT(  T  )/PW 

3170=  WRITE!*, *)'!=', I,'  PERIOD!! >=', PERIOD! I ) 

.3180=75  CONTINUE 

>190=30  READ! 10, *> START 

3200=  DO  3S  Jl, 50 

5210-  RA i 10! 1 )=C (START ) 

5  73(1=  S  I AR  I  =STAR  1  +1 

3.' 40"  5  j  i  UNTINUT' 

3250=  Rt AD! 10,*)  I 

3260-  IF! J.EQ.OUO  TO  60 

3270=  WRITE!*,*) 'ELEMENTS  OF  PHASE  CODE  ARE  GAUSSIAN  SHAPE) 

3280=  MKAN=PW/!2*NE1.EM) 

.5290=  STDEV=MF AN/3.0 

3300=  VAR=STDF.V*STDEV 

.5310=  GAUSS=  1 . 0/SQRT !  2 . 0*  VAR*3 .1416) 

3330=  GFACT=0 .592 358 1  * ! PW/NF.I  Lrt ) *( PW/NEt.EM ) 


3330-  WR I TF  <  * , 305 ) MEAN , STOEV , VAR , GF AC  T 

3340=  TAL  LY=0 

3350=  00  55  IA=1,NPRT 

3360=  CGIJNT= TALLYfl 

3370=  00  45  1=1 ,NFLEM 

3380=  HO  40  J=COUNT ,COUNT+I TIME-1 

3390=  T I ME  =  <  PW/ 1 T I ME ) * ( J-COUNT + 1 ) /NELEM 

3400=  POWER=T I ME-MEAN 

3410=  WT=GAUSS*EXP<-P0WFR#P0WFR/(?.0*VAK) ) 

3410=  AM  1 2=CMPLX (0.0, RAT  10 ( I ) *3 . 1 416  > 

3430=  1 1  .j)=WT*CEXP<AMT?> 

3440=40  CONTINUE 

3450=  COUNT=COUNT+IT IMF 

3460=45  CONTINUE 

3470=  W)  50  J=COUNT .COUNT FPFR]OD(  I A) -SUM  1 

3480=  2<J)=<0. 0,0.0) 

3490=50  CONTINUE 

3500=  T At  I.  Y=TALLY  FPFRIOLK  In  > 

3510=55  CONTINUE 

35.0=  GO  TO  100 

3530=60  WRITE**,#) 'FLEMFNTS  OF  PHASE  COHE  ARF  RECTANGULAR 

3540=  ('SHAPED ' 

3550=  GFACT=1 .0 

3560=  TALI.Y=0 

3570=  AST ART=1 

3580=  l\FNO-NFI  FM 

3590=  DO  80  IA=1 »NPRT 

3600=  COUNT3 TAL  I  Y+l 

3610=  DO  70  I=KSTART,REND 

3620=-  DO  65  J=COUNT, COUNT HTTME-1 

3630=  AM  T 2=CMPLX ( 0 , 0 ,RAT ID  U ) *3 . 1 4 16 ) 

3640  •  /< J)=R<J)*CEXP(AMT2) 

3650=65  CONTINUE 

3660=  COUN  T =COLINT+ 1 T I ME 

3670=70  CONTINUE 

3680=  DO  75  J=COUNT,COUNT+PERIOn( JA)-SUM-1 

3690=  ZCJ)  =  *0. 0,0.0) 

3700=75  CONTINUE 

3/10=  KSTART=KENIiEl 

3720=  t\END=KEND+NELEM 

3730=  1 F ( KEND • GT . 50 ) THEN 

3740=  KS(ART=1 

3750=  KFNO=NELEM 

3760=  END  IF 

3770=  TALI  Y=TAU.Y+PERIOD< IA) 

3780=80  CONHNUE 

3790=  mi  85  1=1, SAMPLE 

3800=  EKFQ(I)=0.0 

3810=85  CONTINUE 

7H?0=  n  SUM#  fOTAI  /PW 

3830=  N  SUM 

3840=100  RETURN 

<050=300  FORMAT-  IX.  FOR  I  =  "  ,  L  .3 , 5X . '  PHASE  =  '  .12,  '  *P  1  - 
3860-305  FORMAT •'  I  X «  ' M- '  ,r'1?.lu.3X,  ST<=  ,F12.10. 


3870=  @!X,  On  fFl2. 10, 3X.  'I3FAL='  ,F12. 10  * 

3880=  E  Nil 

Sb90- 

3900=**** *** t************************************* 
3910=*****+* M ****** ****************************** 

3920=****  3UBR0UTINE  FM 

3930=****  WRITTEN  BY  KENNETH  W.  ALBERT/ i REF  IRKED; 

3940=****  i8  OCT  83 

3950 -**********>************** »*************»>**** 
3960=********» ******* ***************************** 

•970=  SUBROUTINE  FH!FREQ,PE£RIOD,NPRT ,PW,N,C,H, lOIAL .PRT > 

3980- 

3990=  INIEOEEi  PER10BU074*, SUM, COUNT, REMAIN, 

•1000=  8 X 1 ,  XI! ,  X3 ,  X4 ,  X5 ,  X6 , ELUDE  ( 58 ) , C  ( 64 ) 

4010-  REAL  FREQ! 1024) ,TxFREQ,PRT!50> 

1020= 

4*».J0- 

4040=  00  1  JJ=1 , 1024 

4050=  rR.£Q<  JJ)=0.0 

4060=1  CONTINUE 

4070=  DF=7.8125 

4080=  00  10  1=1,58 

4090=  Xl=C<I) 

4100=  x?=cmi> 

4110=  X3C!IE2) 

4 J  20=  A4=r<I+3) 

4130=  X5=C!TE4> 

4140=  X6=C (1+5) 

4150-  IF iX?»EO» 1 )X2=2 

4 J  o0=  IF  < X3.F0. 1 )X3=4 

4170-  IF(X4,F.«,1>X4*8 

4180=  IK !X5,KU, 1 )X5=16 

4190=  If  •,X6.EQ.1)X6=32 

4200=  IF ( ■ XI +X24X3+X4+X5+XA  > ♦ EO . 0 ) THE  N 

4710  WRIT!  <*,*>  'ALL  ZERO  CODEWORD  FOUND  IN  SUBROUTINE 

4220=  8  EM, FCODE!*)  SET  FQUAl  TO  64' 

4230=  FCQDE< I )=64 

4240=  END  IF 

4250=  FCODE! I )=XtEX2fX3EX4EX5EX6 

4260=10  CONTINUE 

4270-  URITF<*,*)  'FCODE  U>  VALUES' 

4280=  DO  7  KR-1,58 

4790=  WRITE!*,*)  FCODE (KR) 

4300=7  CONTINUE 

4310=  READ  ( 10,*) ISUBEl 

4320=  II IMK=N/ISUBEI 

43  30=  M- llIrtF*ISUBEL 

4340=  DO  5  J=1 ,NPR  T 

4350=  PERIOD! J)=N*PRT(J)/PW 

4360=  WRITE!*,  *)'!  =  ',  J, '  PERIOD!  1  )  =  ',  fl  K  Hilt!  .1 ) 

4370=5  CONTINUE 

4380=  rt=N*T(ITAl  /PW 

4390=  8AND=500 

4400=  WRITE!*,*) 'BANDWIDTH  SWEPT  =', 


:  17 


BAND, 'HZ 


4 


4410=  RF.AD<10,*)K 

44?0-  IFlK.EQ.ll  60  TO  25 

4430=  IFIR.EQ.2)  60  TO  45 

4440=  THK.EQ.3)  60  TO  60 

4450=  IF(K.EQ.4)  60  TO  75 

4460=  [F(K.FQ.S)  60  TO  VO 

4470=  WRITE!*,*) 'LFN  IMPt  EMENTED' 

4480=  sun=o 

44V0=  HO  20  1=1 »NPRT 

4500=  CQUNT=SUM+.1 

4510=  REMAIN=COUNT +N- 1 

4520=  DO  15  J=CQUNT, REMAIN 

4530=  TIME=(PW/N)*( J-SUH-1 > 

4540=  FREQ<  J ) =BAND*T IMfc  *  T TME*3 . J  4 1 6/PU 

4550=15  CONTINUE 

4560=  SUM=SUM+PFKIOD(I> 

4570=20  CONTINUE 

4580=  60  TO  200 

4590=25  WRirE<*.*>'VFM  IMPLEMENTED' 

4600=  SUM=0 

4610=  DO  40  1=1 ,NPRT 

4620=  COUNT =SUH+t 

4630=  RFnAlN=COUNTFN-l 

4640=  DO  30  J=COUNT, REMAIN/;' 

4650-;  TIMF-PW/N*<N-J*2-SUM> 

4660=  FREQ(  I)  -0,5»BAND*TIME*T]ME*3. 141 6/PW 

4670=30  CONTINUE 

4660=  DO  35  J= REMMIN/2+i .KFnMIN 

4690=  TIME=PW/N*< J*2-N-SUM-1 ) 

4700=  FRFCK J)  0.5*BAND*TIMfc*TIME*3, 1416/FW 

4710=35  CONTINUE 

4/20=  SUh=SUM+PFRIOD(  I ) 

4730=40  CONTINUE 

47  40  60  10  200 

4/ 50=45  Ui-::li-Uf#)'  INTERPULSE  FH  IMPLEMENTED 

4  760=  3UM*0 

4770=  DO  55  T-l ,NPRT 

4/110-  rnuNr--2,UM+l 

4-90=  hi  MAIN  COUNT +N-1 

4H0»»=  DU  50  J-t .HUNT, REMAIN 

4610=  T  I  ME-- (PW/N  )*<J- 1 ) 

48?0=  FRFQ (  J )  =FCODF  (  D  *DHt  I  ThF  /  ( ?*\j .  1 4 1 6 1 

4(4  40-50  (UNTTNIJF 

4840=  SUM "SUMTPFH (0D< f  > 

4850=55  CONTINUE 

4860=  00  TO  200 

48/0=60  UR (IE (*,*)' TNT RAPUI  SE  FH  TMPLt  rON  TC 1* 

4880=  SUN=0 

4890=  DO  70  .1  =  1  .NPRT 

4900=  C0UNT=SUM+1 

4910=  REHAI N=CQUNT +N- 1 

4920=  NSET=1 

4930=  ISFT-0 

4940=  DO  65  J=LOUNT .REMAIN 


*  r  yv 


4V50=  TIHE=<PU/N)*<J-1> 

4960=  IF  < ISET.EQ, ITIME)THEN 

4970*  KS£T=KSET41 

4980=  IS£T=0 

4990=  END  IF 

5000=  IS£T=ISET+1 

5010=  £R£Q!J>=FC0DE!KSET)*DF*TIH£/!2*3.1416> 

5020=65  CONTINUE 

5030=  SUH=SUH*P£RIOD! I ) 

5040=70  CONI INUE 

5050=  60  TO  200 

5060=75  UNITE!*,*) ' INTERI  UI.5E  FH  WITH  INTRAPUL3F  CHIRP 

5070=  SUH=0 

5000=  DO  85  1=1 ,NPRV 

5090=  COUNT=SUMH 

5100=  RErtrtIN=CnUNF+N-l 

51 10=  DO  80  J=COUNT, REMAIN 

5120=  T I  ME = ! PW/N )  *  !  J-  SI  IM- 1  > 

5150=  FREQ! J)=!FCOO£! 1 )/64) XBANDXTIMEXTIMEX3* 1416/PW 

5140=80  CONTINUE 

5150=  SUM=SUM+PERIOD< I ) 

5160=85  CONTINUE 

5170=  GO  TO  200 

5180=90  WRITE!*,*)  'INTRAf  tll.SE  FH  U/SUBPULSE  CHIRP' 

5190=  SUA=0 

5200=  DO  TOO  I=l,NPRr 

5210=  CQUNT«SUHH 

5220=  R'fc.MA  f  N=CUUN  f  FN-1 

5230=  Rser=i 

5240=  ll>0 

5250=  \SET=0 

5260=  HU  95  j=COUNI, REMAIN 

5270=  fiHE  =  !PU/N)*!  J-S'JM-IC-i) 

5280=  IF  !  <  J-SUH-IC)  .Ed.  ITIMF )  IOiC+ITIME 

5290=  l  F  •.  I  SET .  EQ .  I T I  HE )  THEN 

5300=  KSET=KSET+1 

5310=  ISET=0 

5320=  END  IF 

5330=  .  lSET=ISETfl 

5340=  FREQ!  J)  =  (FCODE  (KSt D/64 ) *BAND*TIMEXT IME*3.  i4i6«‘PW 

5350=95  CONTINUE 

5360=  SUM=SUM+PERIOIi!I) 

5370=100  CONTINUE 

5380=20!;  RETURN 

5390=  END 

5400=***********#***********************¥*****¥**************** 
54\0=***¥********X*******¥***#********************************.«t 
5420=****  SUBROUTINE  RfcXPL 

5430=****  WRITTEN  BY  KENNETH  W.  ALBERT / ( REF l REED ) 

5440=****  1 8  OCT  83 

5450=*******************************************************m 
5460=***¥*¥*******¥¥*¥****< *********¥¥**¥*¥*¥**¥*¥*¥*<*+ ******* 
5470=  SUBROUTINE  RECPL<PER10D,Z,GFACT,NPRr ,N) 


5490-  COMPLEX  Z(1024) 

5500=  INTEGER  SUM, COUNT, RENA I N,PFRT OIK  1024) ,F 

5510=  DOUBLE  PRECISION  DSEED 

5520= 

5530=  SUM=0 

5540=  GFACT=1.0 

5550=  WRITE !*,*) 'PULSF(S)  IS  RECTANGULAR  SHAPFU' 

5560=  READ!10,*)F 

5570=  IF(F ,FQ.O)  GO  TO  20 

5580=  WRITE!*,*; 'NQN-FLUCTUAT iNb  TARGET 

5590=  DO  15  1=1 ,NPRT 

5600=  C0UNT=SUHF1 

5610=  KF  rtA JN=f GUNT  +N  - 1 

5620=  DO  5  J=COUN  ? , RENA  IN 

5630=  Zi  L*--(  1 .0,1.0) 

5640=5  CONTINUE 

5650=  DO  10  J=REMA.TN+J ,  SUlitRERlGO*,  1  > 

5660^  7(J)-M0, 0,0.0) 

5670=10  CONTINUE 

5680=  SUM=SUH+PERIUO( I ) 

5690=15  CONTINUE 

5/00  -  «>!)  TO  40 

5710=20  WRITE!*,*)  RAVI  E  T OH  AMPLITUDE  FI  UCJUATIONV, ' 

5715=  l)SFFD=69.D0 

5770=  DO  35  I- l ,NPRT 

5770=  CnMNr=SUHEl 

5780=  RENA 1 N=COUN  T +N- 1 

5790=  00  25  J=C0UNT.RE»4AIN 

5/V2  =  rt  GUIDES!  LIKE  ID) 

5794=  B=  *ALOG! A) 

5796=  COF=SQRT!B) 

5800=  z(.i>=coe 

5810=75  CONTINUE 

5870  -  DU  30  J=RFMAIN+1 ,  SIJMEPERIUD!  I ) 

5830=  Z<J)--(0. 0,0,0) 

5840=30  CONTINUE 

5850=  SI  JM=SUMf  PERIOD!  I ) 

5860=35  CONTINUE 

5870=40  RETURN 

5880-  END 

5890=********************************************»* ********* 
5900=*****  ************************************  4  ************* 
5910=*****  SUBROUTINE  SETUP 

5V.’0=*****  WRITTEN  BY  KENNETH  W,  ALBER I  /  (REF  I  REED ) 

5930=*#***  18  OCT  83 

j'/40=##*  ***##*#*###*'#*##***#*#**#**'***#**##*##**#**  ********* 
5950=****###* I*######*####*##*###*######*#######* *********** 

5960=  SUBROUTINE  S£TUP(D,Z,FREQ,XVALS,YVALS,UALUL,TOTAl  .SAMPLE, M, 

5970=  @1WK,WK,PW,N) 

5980  ; 

5990=  COMPLEX  XVA1..S ( 1024 ) ,  YVALS ( 1024 ) , Z ( 1 024 ) , UAL UF , AM 1 1 . MM  T2 . 

6000=  8C TIME, AMT3, FILTER! 1074  > ,AHT4 

6010=  INTEGER  SAMPLE.NSCAT  ,1) 

6020=  REAL  FREQ! 1074) 


6030= 
6040= 
6050= 
6060 = 
6070= 
6080= 
6090= 
6100= 
6110= 
6120= 
6130= 
6j40= 
6150= 
6160= 
6170= 
6100= 
6190= 
6200= 
6210= 
6220= 
6230= 
6240= 
6250= 
6260= 
6270= 
6280= 
6290= 
ft  300= 
6310 
c,320= 
6330= 
6340= 
6350= 
6360= 
6320= 
6300= 
6390= 
6400= 
6410= 
6420= 
6430= 
6440= 
6450= 
ft46<>= 
6470= 
ft480  = 
.-.490- 
6500- 
6510 
.•>520 
6530 
6540 


on  5  I=M+1 .SAHPl F 
TVALS< l)=CMPLX(0, 0.0.0) 

5  CONTINUE 

DO  10  I=Nf  1 .SAMPLE 
XVALS(I)=CMPIX<0.0,0.0) 

10  CONTINUE 

DO  15  1=1, M 

AMT2=CHPLX (0.0. FREQ ( I ) > 

TIME=TQTAL/M*<I-1> 

CT  IME=CMPLX<TIME,0.0) 

AhT3=VAl  Ut*CTIMEtAHT  2 
YVALS ( 1 ) =Z ( I ) *CEXP ( AMT3 > 

15  CONTINUE 

DO  20  1=1, M 
AMT 1=Z( 1+(M-I ) ) 

AM  T?=CMPLX  < 0.0,  -FREfi  < 1  +  <  M-I >  >  > 

XVAl S<I)=AMT1*CEXF<AMT2) 

20  CONTINUE 

IF(D.EQ.O)  GO  TO  200 
READ< 10,*) ALPHA 
DO  40  J=l, SAMPLE 
TIrtE=(PU/N)*( J-l) 

AMT4=CrtPLX( -ALPHA4T  Trtt  *T IrtE ,0.0) 

FILTER! J)=CEXP(AHT4  > 

40  CONTINUE 

CALL  FFTCC ( YVAL S , SAMPLE , 1 WK , UK ) 

DO  50  1=1, SAMPLE 

YVrtl S ( l ) = YVALS < l ) *F ILTER ( I ) 

50  CONTINUE 

DO  51  1=1, SAMPLE 

YV6I  S  ( I )  =C(IN.JG ( YVALS ( I ) ) 

51  CONTINUE 

CAL L  FFTCC ( YVALS, SAMPl  E , I WK , WK ) 

DO  52  1=1, SAMPLE 

YVALS ( I )=CONJG< YVALS ( I ) ) /SAMPLE 

52  CONTINUE 

200  RETURN 

END 

************************************************** 
************************************************** 
****  SUBROUTINE  OUT 

****  WRITTEN  BY  KENNETH  W.  ALBERT / < REF i REED ) 

**#*  1R  OCT  83 

y************************************************* 
****«****«******¥«********************»*********** 
SUHRI  HIT  INF  OUT  <  MAG , NUM , PW , TOT AL , h , D ) 


RFAI  MAG (1024) 
INTEGER  D 


UR  j  1E(*.3G0)PW,NIIM 
IMO.EU.O)  THEN 
WHITE  (*,.305) 


•  .VV'.-'' 
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